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ABSTRACT 


Fortran  computer  subroutines  have  been  developed  for  machine 
plotting  of  radar  or  radio  vertical-plane  coverage  diagrams  for  the 
case  of  interference  between  direct  and  sea-reflected  waves.  One  set 
of  subroutines  plots  detection  or  constant -signal-level  contours  on  a 
range-height-angle  chart.  A  second  set  plots  signal  level  in  decibels 
as  a  function  of  range  for  a  target  at  a  fixed  (low)  altitude.  The  first 
type  of  plot  is  valid  for  antenna  heights  that  are  within  a  few  hundred 
feet  of  the  water  and  for  targets  that  are  at  much  higher  altitudes  and 
not  too  close  to  the  horizon.  The  second  type  is  valid  for  antenna  and 
target  heights  that  are  both  less  than  a  few  thousand  feet  altitude,  with 
no  restriction  on  minimum  altitude.  Normal  atmospheric  refraction  is 
taken  into  account,  but  scattering  and  ducting  are  assumed  to  be  of 
negligible  effect.  The  frequency  range  considered  is  from  about  100 
MHz  to  10  GHz.  The  factors  taken  into  account  are  frequency,  antenna 
and  target  heights,  antenna  beamwidth,  tilt  of  the  beam  with  respect  to 
the  horizon,  roughness  of  the  sea  (wave  height),  polarization  of  the 
radio  waves,  and  the  calculated  or  assumed  free-space  range.  Sample 
plots  are  shown  and  discussed.  The  plotting  techniques  employed  are 
briefly  described,  and  listings  of  the  essential  Fortran  subroutines  are 
given  in  an  appendix. 


PROBLEM  STATUS 

This  is  an  interim  report  on  a  continuing  NHL  Problem. 


AUTHORIZATION 

NRL  PROBLEM  R02-64 
Project  SF-11- 141 -005-6173 


Manuscript  submitted  March  16,  1970. 


ii 


8 


MACHINE  PLOTTING  OF  RADIO/RADAR 
VERTICAL-PLANE  COVERAGE  DIAGRAMS 


INTRODUCTION 

When  the  antenna  of  a  radar  or  radio  station  in  the  VHF  to  microwave  frequency 
range  overlooks  a  reflecting  surface,  such  as  the  sea,  the  free-space  radiation  and  re¬ 
ception  patterns  are  modified  t>y  the  presence  of  a  reflected  wave,  in  addition  to  the  wave 
that  goes  by  a  direct  path  from  the  transmitter  to  the  receiver  (by  way  of  the  target  in 
the  radar  case).  The  modified  pattern  in  the  vertical  plane  is  the  result  of  phasor-vector 
addition  of  the  direct  and  reflected  waves;  it  is  an  interference  pattern,  withlalternate 
maxima  and  minima,  or  lobes  and  nulls.  In  optics  this  phenomenon  is  known  as  the 
Lloyd's  mirror  effect.  It  is  an  effect  well  known  to  naval  radar  operators,  and  the  math¬ 
ematics  of  the  problem  have  long  been  well  understood.  In  the  lobe  maxima,  the  radar 
echo  signal  may  be  increased  by  as  much  as  12  dB  compared  to  the  free-space  signal. 

In  the  minima,  or  nulls,  on  the  other  hand,  the  signal  strength  may  theoretically  go  to 
zero. 

It  is  frequently  desired  to  make  a  plot  of  the  vertical -plane  coverage  diagram  of  a 
radar,  taking  into  account  this  interference  effect,  in  order  to  predict  regions  in  which  ; 
targets  will  and  will  not  be  detected.  The  computations  Required  for  such  a  plot,  and  the 
actual  plotting,  are  extremely  tedious  even  though  basically  straightforward.  The  obvi¬ 
ous  present-day  solution  of  this  type  of  problem  is  the  use  of  digital  machine  computation 
and  machine  plotting.  The  purpose  of  this  report  is  to  describe  a  procedure  of  this  type 
that  has  been  developed,  and  to  show  some  typical  resulting  plots.  It  is  believed  that  this 
procedure  constitutes  a  valuable  tool  for  the  system  designer  and  operational  analyst.  It 
allows  detailed  plots  to  be  produced  in  a  short  time  and  at  moderate  cost.  These  plots 
enable  the  effects  of  various  parameters  to  be  studied  visually.  To  attempt  the  same 
thing  by  manual  methods  would  take  an  inordinately  long  time  and  would  be  prohibitively 
expensive. 


INTERFERENCE  OF  DIRECT  AND  REFLECTED  WAVES 
Pattern-Propagation  Factor 

The  geometry  of  the  interference  problem  is  shown  in  Fig.  1,  with  heights  greatly 
exaggerated  relative  to  range.  The  mathematical  formulation  of  the  problem  that  follows 
is  a  modification  and  slight  generalization  of  the  one  given  by  Kerr  (1),  except  for  Eqs. 
(24),  (25),  and  (26),  which  are  derived  in  Appendix  A. 

It  is  apparent  that  the  direct-path  wave  and  the  reflected  wave  traverse  different 
distances  in  going  from  the  antenna  to  the  target;  their  path  difference  is 

8  -  R,  +  Rj  -  R  ■  (1) 

The  phase  difference  of  the  two  interfering  waves  due  to  this  path  difference  is  2 vb/k  ra¬ 
dians,  where  \  is  the  radar  wavelength.  The  total  phase  difference  a  of  the  two  wayes 
is  this  path-length  phase  difference  plus  the  phase  change  4>  that  occurs  in  the  reflection 
process. 
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Fig.  1- Geometry  of  the  surfaoe- reflection 
interference  problem 


If  the  two  interfering  waves  have  (at  least  approximately)  the  same  vector  (spatial) 
direction,  then  the  most  important  factor  in  thoir  addition  is  their  phasor  relationship. 
The  resulting  electric -field  phasor  is 


E  ,  |Ed|  ♦  I  Er|  e",a  .•  (2) 

where  Ed  is  the  electric  field  of  the  direct-path  wave,  Er  is  that  of  the  reflected  wave, 
and  a  is  their  phase  difference. 

Kerr  (1)  defines  the  propagation  factor  F  as  the  ratio  of  the  actual  field  |E|  at  a 
point  in  space  to  that  which  would  exist  in  free  space  at  the  same  distance,  |  E0 1 .  If  |  Ed  | 
is  assumed  equal  to  |  E0 1 ,  as  it  will  be  in  the  absence  of  absorptiomlbsses  or  abnormal 
refractive  effects,  from  Eq.  (2)  it  is  apparent  that 


If  the  path  difference  5  is  assumed  to  be  very  small  compared  to  R  (as  will  always 
be  true  in  situations  of  practical  interest),  for  a  flat  smooth  surface  the  ratio  |Er|/|E0| 
is  equal  to  the  magnitude  p0  of  the  reflection  coefficient  of  the  reflecting  surface. 
Therefore 


F  =  i  1  *  P0  e“Ja|  =  1/1  +  p2  +  2p0  cos  a  |  .  (4) 

If  the  surface  is  not  smooth  and  flat,  p0  is  replaced  by  p' ,  which  is  the  product  of 
three  separable  factors:  a  factor  dependent  on  the  roughness  of  the  surface,  one  that  de¬ 
pends  on  its  curvature  (e.g.,  jphericity  of  the  earth),  arid  one  that  expresses  the  inherent 
reflectivity  of  the  material.  If  the  reflection  coefficient  is  p0  for  a  plane  smooth  sur¬ 
face,  a  roughness  factor  r  and  a  divergence  factor  D  can  be  defined  to  take  into  account, 
respectively,  the  roughness  and  the  curvature.  Then 

p‘  -  r lipQ  .  (5) 

Thus  far  F  has  been  defined  as  simply  the  "propagation  factor."  It  is  necessary, 
however,  to  take  into  account  in  this  factor  that  the  antenna  vertical-plane  pattern  can 
also  cause  a  difference  in  strength  of  the  direct  and  reflected  rays.  The  factor  f  (&)  de¬ 
fines  the  antenna  pattern.  It  is  the  ratio  of  the  field  strength  in  direction  e  relative  to 
that  in  the  beam-maximum  direction,  where  e  is  angle  measured  with  respect  to  an 
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arbitrary  reference  direction  in  the  vertical  plane.  If  (Fig.  1)  is  the  direction  of  the 
direct  ray  as  it  leaves  the  antenna  and  is  the  direction  of  the  reflected  ray  at  the 
same  point,  it  may  or  may  not  be  true  that  f  (8t)  -  f  (d2).  It  is  therefore  necessary  to 
define  the  pattern  -propagation  factor  (which  F  will  henceforth  ienote)  as 


F  *  |f(5,)  +  e“J“| 


f<*,) 


1  t  p‘ 


f  (0a) 
f  <**> 


(6) 


As  thus  defined,  F  is  the  ratio  of  the  actual  field  strength  at  a  point  in  space  to  that  which 
would  exist  at  the  same  distance  from  the  antenna  in  free  space  in  the  beam  maximum. 

A  generalized  reflection  coefficient,  x,  can  now  be  defined: 


,  ft**)  rDPof(02) 
X  =  P  f<5,)  r  f(0t) 

The  pattern  propagation  factor  then  becomes  simply 

F  =  f (&t)  jl  +  xe'ja[ 


(7) 


in  which 


-  f(0,)  j/l  +  "xT  +  2x'  cos  a | 
2vb 

a  -  -  f  4>  , 


(8) 

(9) 


where  <t>  is  the  phase  angle  of  the  complex  reflection  coefficient.  This  phase  angle  is 
defined  as  the  phase  lag  of  the  reflected  wave  relative  to  the  phase  of  the  incident  wave. 


Path  Difference  and  Divergence  Factor 


If  the  earth  were  flat,  it  is  readily  shown  that  in  the  terms  of  Fig.  1, 


s 


2hih» 

d 


(10) 


where  d  is  the  horizontal  distance  between  the  antenna  and  the  target.  This  result  is 
obtained  by  assuming  that  dJ  »  hj3 .  The  slant  range,  r,  is  given  by 


R  =  V(hj-hj)1  +  d3  .  (11) 

When  the  heights  h,  and  ha  and  the  distance  d  are  such  that  the  earth's  curvature 
is  significant,  Eqs.  (10)  and  (11)  are  not  really  correct,  although  in  many  practical  situa¬ 
tions  they  are  excellent  approximations.  The  exact  curved-earth  solution  for  b  is  com¬ 
plicated,  because  it  requires  solution  of  a  cubic  equation.  The  method  that  has  been  used 
for  the  computer  programs  being  described  is  that  given  by  Kerr  (1),  p.  113  ff,  in  terms 
of  a  pair  of  parameters  s  and  T.  These  parameters  are  defined  by  the  following 
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expressions,  in  which  h,  and  h2  axe  now  (and  henceforth)  defined  as  heights  above  the 
earth's  surface  rather  than  above  the  tangent  plane  as  in  Fig.  1: 

s  -  - — - - . 

■  \/  2  ap  ( •J h  j  +  \/hj ) 

T  i/h j or  \J h 2 /h ,  , 

in  which  ae  is  the  "effective  earth's  radius,"  for  the  assumed  conditions  of  atmospheric 
refraction,  and  the  choice  of  definition  of  T  is  made  to  result  in  T  <  1 .  The  effective 
earth's  radius  concept  is  based  on  the  treatment  of  atmospheric  refraction  given  by 
Schelleng,  Burrows,  and  Ferrell  (2);  it  is  an  approximate  method  which  is  adequate  for 
the  situation  being  considered,  for  antenna  and  target  heights  not  in  excess  of  a  few  thou¬ 
sand  feet.  The  accepted  value  of  ne  for  "standard  refraction"  is  4/?  times  the  actual 
earth's  radius.  If  the  actual  radius  is  taken  to  be  3440  nautical  miles,  and  if  d  is  in 
nautical  miles  and  h,  and  h2  are  in  feet,  Eq.  (12)  becomes: 

s 


In  terms  of  these  parameters, 
the  divergence  factor  D  of  Eq.  (5). 
formula  for  5  so  that  it  becomes 


Kerr  presents  curves  for  J  and  D  as  functions  of  s  and  r.  To  compute  these  quan¬ 
tities,  it  is  necessary  to  employ  the  intermediate  parameters  s,  and  S2,  which  are  re¬ 
lated  to  S  and  T  by 


1.2287  (v/h7  t  yTij) 


(14) 


a  correction  factor  J  can  be  computed,  as  can  also 
This  correction  factor  is  applied  to  the  flat -earth 


2h.h2 

s  * - -  J  (S.T) 


(15) 


(12) 

(13) 


S  -•  (S,  *  S2T)  (1*T)  .  (16) 

Sj  j/(l  S,2)J  i  4  S j2  T2  i  S,2  ■  lj/(2S.T)  .  (17) 

An  iterative  procedure  is  used  in  machine  computation  to  find  values  of  s,  and  s2 
that  simultaneously  satisfy  Eqs.  (16)  and  (17)  for  given  values  of  s  and  T.  When  s,  and 
s2  have  thus  been  found,  j  and  D  are  given  by 

J  ( 1  -  Sj2 )( i  -  s22 )  ,  (18) 


n 


4S,2  S2T 

sd:V>(i  *  ‘O 


(19) 


The  iteration  is  performed  using  tire  computer  subroutine  named  INVERT.  A  listing 
of  the  Fortran  statements  comprising  this  subroutine  is  given  in  Appendix  B. 
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Surface-Roughness  Factor 

The  roughness  factor  r  of  Eq.  (5)  is  computed  as  a  function  of  the  sea- wave  height 
H,  the  radar  wavelength  K,  and  the  grazing  angle  4>,  Fig.  1,  using  a  formula  given  by 
Ament  (3):* 


r  -  ?xp 


-2 


( 


2v  H  sin 
\ 


(20) 


Beard,  Katz,  and  Spetner  (4)  have  conducted  experiments  whose  results  are  in  approxi¬ 
mate  agreement  with  this  formula.  It  was  assumed  in  deriving  Eq.  (20)  that  the  sea  sur¬ 
face  has  a  Gaussian  height  distribution  and  that  H  is  the  standard  deviation.  For  sea 
waves  that  have  approximately  sinusoidal  shape  H  can  be  equated  to  H  7(2/2),  where  H' 
is  the  crest-to-trough  wave  height,  since  the  standard  deviation  of  a  sine  wave  of  zero 
mean  value  is  given  by  its  rms  value,  which  is  1/n/T  times  the  amplitude  of  the  wave  or 
1/(2  7*2)  times  the  crest-to-trough  height. 

The  grazing  angle  4>  is  given  by  the  formula 

tan  <p  -  - — -j  K(S,T)  .  .  (21) 


where  the  first  term  is  the  flat-earth  formula  and  K  is  a  curved-earth  correction  factor 
given  in  terms  of  S„  s2,  and  t  by 


% 


K  -- 


(1  -  Sj2)  +  T2(  1  -  S2) 
1  i  T2 


(22) 


Approximate  Method  for  Distant  Targets 

The  spherical-earth  computation  of  s  and  V'  using  the  parameters  s  and  T,  as  just 
described,  requires  a  priori  knowledge  of  the  target  height.  For  plotting  coverage  dia¬ 
grams  it  is  convenient  to  assume  a  target  elevation  angle  rather  than  a  height.  The  1 
pattern-propagation  factor  can  be  computed  as  a  function  of  angle  if  the  assumption  is 
made  that  the  target  is  at  a  range  much  greater  than  the  distance  from  the  antenna  to  the 
reflection  point.  The  method  of  doing  this  for  the  flat-earth  assumption  is  well  known; 
the  path  difference  in  tuat  case  becomes  simply 

S  -  2h,  sin  0A  .  (23) 

This  result  can  be  deduced  from  Fig.  1  by  assuming  that  the  earth  is  flat  and  the  target 
is  so  distant  that  <p  =  0d.  (The  lines  labeled  R  and  R,  are  then  considered  to  be  parallel.) 

It  is  possible  to  extend  this  method  to  apply  to  tire  spherical-earth  case.  The  result¬ 
ing  equation  for  &  is 


5  -  i/i//  *  «,(ne  t  h)  2  sin3  (#d  *  <p) 


(24) 


*Ament  states  that  this  formula  was  derived,  Independently,  by  Rekerts  and  by  MacFarlane,  during 
World  War  U,  / 


() 
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wh»re 

>/■  \/(»nn  l]3  1  Jh,  "  ''»>  i^,|)  3  (25) 

unit  la  (he  efiacilvo  earth’s  radius  as  In  Eq.  (12).  Tho  grazing  angle  v  lf>  given  by 

r  '  <\\  i  ■>  ■  (26) 

The  derivation  of  these  results  Is  given  In  Appendix  A. 

The  calculation  of  the  divergence  factor,  n,  ear  also  be  approximated  in  terms  of  <■,, 
by  a  formula  given  by  Kerr  (1),  Eq.  (474),  p.  137: 


r 

/ 

l 

x 

\ 

.iL 

.1 

v^'n) 

where 


f  'C  " 

v  1", . . 


Uofloctlon  Coefficient  of  Sea  Water 

Tho  most  common  situation  that  gives  rise  to  iin  Interference  pattern  is  reflection 
from  the  sea;  therefore  It  Is  of  Importance  to  evaluate  the  intrinsic  reflection  coefficient, 
and  the  roUeotion-coefficlent  phase  angle  ,/■  for  sea  water.  The  equations  as  given 
hy  Kerr  (1),  p.  396  ff,  for  the  horizontal-polarization  coefficient  ih  and  the  vertical- 
polarization  coefficient  l  v ,  are 


-u  sin  V  “  \fer  ~  COS2  ./• 

|>  '  '\><h>  «'  '  - 

sin  V  1  v t  ~  cos'  s/' 


V  '  o  (  V  ) 


f=0  sin  1 1’  -  \f f. t.  -  cos2  \ji 
*:c  sin  v'  4  /«c  -  cos2  7/ 


In  which  ^  Ib  the  grazing  angle  (Fig.  1)  and  cc  is  the  complex  dielectric  constant  of  the 
reflecting  material,  given  by 

1  *c  «  I  j ,  J  <  1  “  j60\rr  .  (31) 

The  real  part  of  ec  is  the  ordinary  dielectric  c  onstant,  k  is  the  wavelength,  and  r.  is  the 
conductivity.  The  factor  60  applies  when  a  is  in  mhos  per  meter  and  \  is  in  meters. 

The  values  of  t ,  and  e2  are  functions  of  frequency. 

Fig-ures  2,  3,  and  4  are  plots  of  <’0(h),  p0(v),  and  <£v  produced  by  machine  compu¬ 
tation  and  machine  plotting,  using  valuec  of  t ,  and  t>  shown  in  the  following  table: 


■wmjw ■  •IH.HWU  V 1 1  "HflWklwWTt. 

4»  ,. 
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f  (MHz) 

ei 

o-  (mhos/meter) 

<1500 

80 

4.3 

3000 

69 

6.5 

10000 

52 

16 

These  values  for  f  §  3000  MHz  are  taken  from  Kerr's  Table  5.1,  p.  398.  However,  the 
10000-MHz  value  of  e ,  is  taken  from  Von  Hippe1  (5).  The  value  at  this  frequency  given 
by  Kerr's  table  resulted  in  plots  that  were  in  poor  agreement  with  similar  plots  given  by 
Kerr  (his  Figs.  5.4,  5.5,  and  5.6).  It  was  therefore  evident  that  these  curves  were  calcu¬ 
lated  for  another  value  of  f ,  at  f  =  10000  MHz.  This  conclusion  is  substantiated  by 
Kerr's  comment  on  page  401:  "The  values  of  e,  and  a  used  for  these  figures  were 
taken  from  Table  5.1,  except  for  3  cm"  (f  =  10000  MHz).  Using  Von  Hippel's  value  of  «, 
at  this  frequency  produced  better  agreement  with  Kerr's  curves.  The  variation  between 
the  values  tabulated  was  approximated  by  linear  expressions,  as  follows: 

«i  ■  80  ) 

L  f  v  1500  MHz  ,  (32a) 

cr  ••  4.3J 

«j  -  80  -  0.00733  (  f  -  1500)  1 

f.  1500  ■■  f  ■  3000  ,  (32b) 

a  -  4.3  t  0.00148  (f  -  1500)  J 

-  <39  -  0.00  243  (f  -  3000)  ] 

}  ,  3000  <  f  <  10000  .  (32c) 

<7  -  6.52  t  0.001314  (f  -  3000)  J 


Fig.  2  -  Magnitude  of  the  sea-water  reflection 
coefficient  for  vertical  polarization,  ) 
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Fig.  3  -  Phase  angle  of  the  sea-water  reflection 
coefficient  for  vertical  polarization,  8V 
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Fig,  4  -  Magnitude  of  the  sea- water  reflection 
coefficient  for  horizontal  polarization,  p0(h , 
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The  subroutine  used  to  do  the  computation  of  p0  find  4>  for  Figs.  2,  3,  and  4  is 
named  SEAREF;  it  is  a  Fortran  coding  of  Eqs.  (29)  through  (32).  A  listing  of  the  sub¬ 
routine  is  given  in  Appendix  B.  Curves  are  not  given  for  4>h>  since  the  variation  from 
4>  =  180  degrees  over  the  range  of  frequency  considered  is  less  than  4  degrees  (<fr,  ~  184 
degrees  at  f  -  10000. MHz  and  ^  =  90  degrees).  Also,  Po(h)  k  1  over  this  range  of  fre¬ 
quency  and  grazing  angle.  (The  variation  is  magnified  in  Fig.  4  by  the  expanded  ordinate 
scale.)  The  subroutine  should  not  be  used  for' frequencies  appreciably  above  10000  MHz, 
because  Eqs.  (32c)  are  probably  not  valid  at  higher  frequencies. 


THE  INTERMEDIATE  REGION 

The  equations  that  have  been  given  permit  computation  of  F  when  there  is  interfer¬ 
ence  between  a  direct  ray  and  a  reflected  ray.  This  situation  exists  when  a  target  is 
above  the  radar  horizon.  The  target  is  then  said  to  be  in  the  interference  region  (Fig.  5). 
Below  the  horizon,  ray  theory  predicts  a  shadow  —  zero  signal  strength.  There  is,  how¬ 
ever,  some  signal  due  to  diffraction;  consequently  the  regions  above  and  below  the  hori¬ 
zon  are  called,  respectively,  the  interference  region  and  the  diffraction  region.  If  the 
earth  had  no  atmosphere,  the  only  signal  in  the  shadow  region  would  be  that  due  to  dif¬ 
fraction.  Actually  there  are  also  contributions  due  to  scattering  by  irregularities  of  the 
high-altitude  atmosphere  and,  on  occasion,  there  can  be  strong  beyond-the-horizon  sig¬ 
nals  due  to  anomalous  refraction  (ducting).  The  plotting  method  described  here  takes 
into  account  only  the  ''normal”  atmosphere,  in  which  there  is  no  ducting,  and  it  does  not 
consider  scattering.* 


Fig.  5  -  Interference,  intermediate, 
and  diffraction  regions 


In  the  absence  of  ducting,  the  field  well  beyond  the  horizon  at  VHF  and  above  is  too 
weak  to  be  of  importance  for  radar.  Detection  can  occur  for  large  targets  that  are  only 
slightly  beyond  the  horizon  with  high-sensitivity  radars.  However,  a  large  increase  of 
power  produces  only  a  very  small  range  increase 4n  this  region. 

In  a  region  very  close  to  the  horizon  but  above  it,  interference  calculations  of  the 
type  that  have  been  described  are  not  valid.  They  are  based  on  the  principles  of  ray 


*The  scatter  mechanism  can  be  the  chief  contributor  to  the  total  signal  at  distances  well  beyond  the 
horizon  (well  below  the  tangent  ray).  At  small  distances  beyond  the  horizon,  however,  the  diffrac¬ 
tion  field  predominates.  At  low  frequencies  and  with  vertical  polarization,  the  surface  wave  can 
also  be  of  significance,  but  it  is  a  negligible  factor  above  100  MHz,  especially  with  horizontal 
polarization. 
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optics,  which  are  not  valid  in  this  region.  A  rigorous  electromagnetic -wave  solution  of 
the  problem  in  this  region  is  very  difficult,  as  discussed  by  Kerr  (1).  The  mathematical 
series  that  describes  the  field  converges  very  slowly.  Further  out,  well  below  the  hori¬ 
zon,  the  series  converges  very  rapidly,  and  the  first  term  adequately  describes  the  dif¬ 
fraction  field.  This  is  termed  a  "one-mode"  solution.  Thus  it  is  possible  to  calculate 
the  field  strength  accurately  (if  anomalous  refraction  and  scattering  are  not  significant) 
in  both  the  interference  and  diffraction  regions  not  too  close  to  the  horizon.  In  the  near 
vicinity  of  the  tangent  ray,  however,  neither  the  interference  calculation  nor  the  one¬ 
mode  diffraction  calculation  are  valid.  Kerr  has  therefore  called  this  the  "intermediate 
region." 

Kerr  proposes  the  following  method  for  estimating  the  field  in  the  intermediate  re¬ 
gion;  he  terms  it  a  method  of  "bold  interpolation."  In  this  method,  a  plot  is  made  of  f 
expressed  in  decibels  (20  log10  F)  as  a  function  of  the  range  R  for  a  target  whose  height 
is  constant,  and  for  a  fixed  height  of  the  radar  antenna.  In  the  interference  region  this 
plot  will  exhibit  the  maxima  and  minima,  or  lobes  and  nulls,  characteristic  of  the  inter¬ 
ference  of  direct  and  reflected  rays.  In  the  diffraction  region,  20  log10  F  decreases 
monotonically,  in  a  quasi- linear  fashion.  Kerr's  method  for  plotting  F  in  the  intermedi¬ 
ate  region  is  simply  to  draw  a  smooth  curve  that  joins  the  curves  validly  calculated  in 
the  interference  region,  sufficiently  above  the  tangent  ray,  and  in  the  diffraction  region 
well  below  the  tangent  ray,  where  the  one-mode  solution  is  permissible. 

The  first-mode  diffraction-field  solution  for  F  can  be  expressed  as  the  product  of 
three  terras  (Kerr,  p.  122): 

F  =  V  (X)  U(Z,)  U(Za)  .  (33) 

in  which  x,  z,,  and  z2  are  the  target  range,  antenna  height,  and  target  height  expressed 
in  "natural  units."  The  natural  unit  of  range,  l,  is  given  by 

L  k,/f  ‘  1  ,  (34) 

in  which  f  is  the  radar  frequency.  For  f  in  megahertz  and  L  in  nautical  miles,  k ,  = 
102.7.  The  natural  unit  of  height,  H,  is  given  by 

H  =  k a/  f 3  3  .  (35) 

For  f  in  megahertz  and  H  in  feet,  k2  =  6988.  (See  Kerr,  pp.  96  and  97,  Eqs.  (351)  and 
(358).)  The  natural  range  x  is  equal  to  R/L,  where  R  is  the  actual  range,  and  the  natural 
height  2  is  h/H,  where  h  is  the  actual  height. 

V(X)  is  called  the  attenuation  factor  and  is  given  by  (Kerr,  p.  122): 

V(X)  =  2v^Te-2°2X  •  (3®) 

or  in  decibels 

VdB  :  10.99  t  10  logJ0  X  -  17.55X  ,  (36a) 

U(Zj)  and  U(Za)  are  called  height -gain  factors,  and  their  calculation  is  very  com¬ 
plicated  (see  Kerr,  pp.  109-112).  However,  a  curve  for  U(Z)  is  given  in  Kerr,  Fig.  2.20, 
p.  128.  In  order  to  be  able  to  use  Kerr's  plot  in  a  computer  program,  empirical  equa¬ 
tions  were  fitted  to  the  curve.  These  equations  give  u<Z)  in  decibels  (udB  •  20  ioBl0  1 j): 
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UdB  =  20  log ,a  Z,  2  <0.6,  (37a) 

UdB=  -4.3  ♦  [51.04  log10  (2/0.6)] 1,4  .  0.6  <  2  <  1  ,  (37b) 

UdB  =  19.85  (2-47-0.9)  ,  2  >  1  .  (37c) 

These  empirical  expressions  fit  the  curves  to  within  about  1  dJB.  These  equations 
have  been  incorporated  into  a  Fortran  subroutine  named  UFCN.  A  listing  of  the  subrou¬ 
tine  is  given  in  Appendix  6. 

These  results  are  valid  for  horizontal  polarization  and  for  the  standard  4/3-earth- 
radius  atmosphere.  Above  100  MHz,  however,  the  differences  in  U(2)  for  horizontal  and 
vertical  polarization  are  negligible  for  z  S  1.  At  low  to  moderate  heights,  this  standard 
atmosphere  is  reasonably  representative  of  the  actual  atmosphere  in  the  absence  of 
ducting. 


LOBE  PLOTTING 

Computer  programs  for  two  kinds  of  plots  have  been  developed,  which  present  the 
interference  lobe  patterns  in  different  ways.  The  first  is  a  plot  of  radar  detection 
contours  in  a  vertical  plane  with  range  and  height  as  coordinates.  The  second  is  a  plot 
of  signal  level  in  decibels  as  a  function  of  range  for  a  target  of  fixed  height. 


Detection-  or  Constant -Signal-Strength-Contour 
Range-Height-Angle  Plots 

The  detection  or  constant-signal  contour  plot  is  made  on  a  coordinate  grid  called  a 
range-height-angle  chart.  The  basic  method  of  computing  this  chart  has  been  described 
(6).  For  the  lobe  plots,  the  chart-plotting  program  has  been  made  a  subroutine,  called 
RHACHT,  such  that  any  maximum  range  and  height  can  be  specified.  The  chart  will  then 
be  drawn  and  properly  labeled  for  these  specified  values.  Then  a  subroutine  called 
LOBES  causes  the  detection  contours  to  be  plotted  to  the  proper  scale  on  this  chart.  The 
method  of  plotting  the  chart  is  such  that  refraction  is  taken  into  account,  on  the  basis  of 
the  CRPL  Reference  Exponential  Atmosphere  with  N,  -  313  (the  surface  value  of  the 
refractive  index  is  0.000313).  (The  chart  could  be  plotted  for  other  values  of  st  by 
changing  two  cards  in  the  RHACHT  subroutine.) 

The  Fortran  parameters  needed  in  the  calls  to  these  subroutines  are: 

XMAX—  the  maximum  x -dimension  of  the  chart,  inches, 

YMAX —the  maximum  y-dimension, 

RMAX  — the  maximum  range  df^the  chart,  nautical  miles, 

HMAX  —  the  maximum  height  of  the  chart,  feet, 

AHFT  —  the  antenna  height,  feet, 

FMHZ—  radar  frequency,  megahertz, 

BWD  —  antenna  half -power  vertical  be  am  width,  degrees, 

WHFT  —  soa  wave  height  (crest-to-trough,  feet), 

TILT  — the  tilt  angle  of  the  antenna  beam  maximum  with  respect  to  the  horizon, 
degrees, 
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THM  —  the  maximum  elevation  angle  to  which  the  plot  is  to  be  made, 

IPOL  —  an  integer,  1  for  vertical  polarization,  2  for  horizontal, 

RFS  —  the  calculated  or  assumed  free-space  range  of  the  radar  on  the  specified 
.  target,  or  for  a  one-way  radio  system,  the  free-space  range  at  which  the 
field  strength  would  have  a  specified  value. 


The  parameters  BVVD  and  TILT  are  needed  for  calculating  the  antenna  pattern  fac¬ 
tors,  f  (0,)  and  f  (02),  Eq.  (6).  The  antenna  vertical-plane  beam  shape  is  assumed,  in 
this  subroutine,  to  be  symmetrical  with  respect  to  the  beam  maximum,  and  to  be  of  the 
"(sin  x)/x"  form,  which  means  that 


f (0)  =  (sin  u)/u  ,  (38) 

where 

u  =  k  sin  0  radians  .  (39) 

The  value  of  k  is  chosen  so  that  f  (0)  a  0.7071  (=  1 /'Ti)  when  6  =  bwd/2.  This 
means 

1  k  =  1.39157/sin  (BWD/2)  .  (40) 

In  the  plotting  subroutine,  account  is  taken  of  the  fact  that  f  (0)  becomes  alternately 
positive  and  negative  for  successive  sidelobes  of  the  antenna  pattern;  when  it  is  negative 
for  the  reflected  ray  and  positive  for  the  direct  ray,  or  vice  versa,  the  effect  is  to  add 
v  radians  to  the  phase  angle  4>  of  the  effective  reflection  coefficient,  x,  Eq.  (7);  that  is, 
x  becomes  negative. 

The  path  difference  s  and  grazing  angle  4>  are  computed  at  the  lower  angles,  where 
earth’s  curvature  is  significant,  by  the  method  of  Eqs.  (24)  to  (27),  and  D  is  computed 
from  Eq.  (28).  The  computation  is  started  at  the  elevation  angle  for  which  the  path  dif¬ 
ference  is  8  =  \/4.  The  path  difference  is  computed  by  both  the  flat -earth  formula,  Eq. 
(23),  and  by  the  curved-earth  formula,  Eq.  (24).  (Since  these  are  both  long-range  approx¬ 
imations,  the  plotted  patterns  are  not  exact  in  the  very -close-range  region,  but  this  fact 
is  not  usually  important.)  When  an  elevation  angle  is  reached  for  which  these  two  results 
differ  by  less  than  0.001  wavelength,  the  flat-earth  formula  is  used  thereafter.  Also, 
the  value  of  D  is  thereafter  no  longer  computed,  but  is  assumed  equal  to  unity,  and  <p  is 
assumed  equal  to  the  target  elevation  angle  0d  (Fig.  1). 

The  basic  method  of  computing  the  radar  free-space  range,  RFS,  is  discussed  by 
Kerr  (1).  The  problem  is  discussed  in  considerably  greatqi*  detail  in  a  recent  two-part 
NRL  report  (7). 

% 

Samples  of  the  plots  resulting  from  this  program,  utilizing  subroutines  RHACHT  and 
LOBES,  are  presented  in  Figs.  6  through  13.  These  samples  show  some  of  the  effects 
that  can  be  studied  by  varying  the  frequency,  antenna  height,  be  am  width,  beam  tilt,  sea 
roughness,  and  polarization,  The  first  plot,  Fig.  6,  for  a  low  frequency,  low  antedpa 
height,  broad  vertical  beamwidth,  smooth  sea,  and  horizontal  polarization,  shows  the 
widely  spaced  well-defined  lobes  that  are  characteristic  of  this  situation,  with  the  lobe 
maxima  reaching  to  virtually  twice  the  free-space  range.  The  second  pl<  Fig.  7,  shows 
the  effect  of  simply  increasing  the  antenna  height,  from  20  feet  to  80  feei  This  causes 
the  lobes  to  be  finer  and  more  closely  spaced,  so  that  there  are  more  of  them,  and  the 
lowest  one  lies  closer  to  the  surface.  The  envelope  of  the  lobe  maxima  is  unchanged. 
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RANGE.  NAUTICAL  MILES 

F  .  100  MHZ  «  ANT  HT  .  20  FT  «  VERT  BH  *  90  DEC  « 

HAVE  H.T  •  0  FT  «  8M  TILT  -  0  OEO  «  FS  RNG  »  100  NM  ■ 
POLARIZATION.  HORIZONTAL 

Fig.  6  -  Detection  contour  plot  for  frequenoy  100  MHz, 
antenna  height  20  ft  above  sea  surfaoe,  vertical  beam- 
width  90  degrees,  smooth  sea,  jero  beam  tilt,  horizontal 
polarization,  and  free-space  range  100  nautical  miles 


F  -  100  MHZ  •  ANT  HI  .  80  FT  .  VERT  BH  ■  90  DEG  ■ 

WAVE  HT  «  0  FT  .  BM  T|LT  «  0  DEC  «  FS  RNC  »  100  NM  . 

POLAR I ZAT 1 0N»  HOR ! ZONTAL 

Fig.  7  -  Detection  contour  plot  for  the  same  parameters 
as  Fig.  6  except  antenna  height  lnoreased  to  80  ft 
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F  »  100  MHZ  «  ANT  HT  =  80  FT  ■  VERT  BU  •  90  DEC  « 

WAVE  HT  -  0  FT  ■  8M  TILT  *  0  DEG  »  FS  RNC  »  100  NM  ■ 

POLARIZATION.  VERTICAL 

Fig.  8  -  Detection  contour  plot  for  the  same  parameters 
as  Fig.  7  except  vertical  polarization 


The  third  plot,  Fig.  8,  is  for  the  same  conditions  as  Fig.  7  except  that  the  polariza¬ 
tion  has  been  changed  from  horizontal  to  vertical.  The  envelope  of  the  maxima  is  now 
radically  changed,  because  of  the  reduced  valufe  of  p0  (see  Fig.  2).  Also,  the  elevation 
angles  of  the  lobe  maxima  are  changed,  with  the  lowest  lobe  at  a  higher  angle;  this  result 
is  due  to  the  difference  in  the  reflection-coefficient  phase  angle  <t>,  which  is  180  degrees 
for  horizontal  polarization  (see  Fig.  3).  The  superiority  of  horizontal  polarization  and 
the  reason  for  its  almost  universal  use  for  VHF-UHF  radar  are  clearly  seen. 

In  the  next  plot,  Fig.  8,  the  polarization  is  again  horizontal,  and  all  other  factors  ex¬ 
cept  the  wave  height  are  the  same  as  in  Fig.  7.  The  effect  of  6-foot  waves  at  100  MHz  is 
seen  to  be  negligible  at  the  lowest  elevation  angles,  but  considerable  above  about  20  de¬ 
grees.  The  result  is  to  reduce  the  magnitude  of  the  reflected  wave  so  that  the  lobe  max¬ 
ima  are  of  reduced  strength  and  the  nulls  do  not  go  to  near -zero  strength,  as  they  do  at 
the  lower  angles  and  for  a  smooth  sea.  At  60-degree  elevation,  whore  the  plot  is  termi¬ 
nated,  only  a  small  interference  lobe  structure  remains;  the  patturn  Is  almost  the  free- 
space  pattern  of  the  antenna,  which  is  3  dB  down  at  45-degree  elevation  angle  (90-degree 
beamwidth),  (The  3-dB  pattern  factor  results  in  a  free-space  range  of  70.71  miles  at 
this  elevation  angle.) 

In  Fig,  10,  the  frequency  has  been  increased  by  a  factor  of  10,  to  1000  MHz,  the 
beamwidth  has  been  reduced  to  10  degrees,  and  the  free-space  range  has  been  decreased 
to  50  nautical  miles.  Other  factors  are  the  same  as  in  Fig.  7.  The  principal  effect  ob¬ 
served  is  the  much  finer  lobe  structure  and  the  much  lower  angle  of  the  lowest  lobe. 

Also  notable  is  the  slightly  reduced  range  of  the  maximum  of  the  lowest  lobe.  This  is 
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Fig.  9  -  Detection  contour  plot  for  the  same  parameters 
as  Fig.  7  except  sea  wave  height  increased  to  6  ft 
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F  »  1000  HHZ  ■  ANT  HT  -  80  FT  ■  VERT  BW  =.  10  DEG  ■ 

HAVE  HT  •  0  FT  ■  BH  TILT  »  0  DEG  ■  FS  RNG  *  50  KIM  ■ 
POLARIZATION.  HORIZONTAL 

Fig.  10  -  Detection  contour  plot  for  the  same  parameters 
as  Fig.  7  except  frequency  Increased  to  1000  MHz,  antenna 
vertioal  beamwldth  reduced  to  10  degrees,  and  free-spaoe 
range  deoreased  to  50  nautical  miles 
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the  result  of  the  divergence  factor,  D,  Eqs.  (5),  (7),  and  (19),  which  is  negligibly  different 
from  unity  for  the  lowest  100-MHz  lobe  but  not  so  for  the  lowest  lobe  at  1000  MHz  (with 
antenna  height  80  feet  in  both  cases).  Also  in  this  plot  the  sidelobe  pattern  of  the  antenna 
is  visible.  The  apparent  increased  fineness  of  the  interference  lobes  in  the  antenna- 
pattern  side  lobes  is  not  real;  it  is  the  result  of  the  nonlinear  angle  scale  of  the  chart, 
which  is  a  necessary  consequence  of  the  fact  that  the  range  and  height  scales  are  differ¬ 
ent. 

The  next  plot,  Fig.  11,  differs  from  Fig.  7  only  in  that  the  wave  height  is  now  2  feet 
instead  of  zero.  Again,  as  at  100  MHz,  this  has  virtually  no  effect  in  the  very  low  angle 
region,  but  it  reduces  the  lobe  strength  at  higher  angles,  and  fills  in  the  nulls,  because  of 
the  reduced  reflection  coefficient.  In  the  first  side  lobe  there  is  only  a  faint  suggestion 
of  a  lobe  pattern,  and  in  the  second  side  lobe  the '’pattern  is  virtually  the  free -space  pat¬ 
tern  of  the  antenna. 

Figure  12  shows  the  effect  of  changing  to  vertical  polarization  at  1000  MHz,  again 
with  a  smooth  sea.  As  at  100  MHz,  the  range  of  detection  in  the  lobe  maxima  is  reduced 
compared  to  horizontal  polarization,  and  the  lobe  angles  are  changed,  but  noticeably  only 
in  the  middle  range  of  angles.  The  very  lowest  lobe  is  only  very  slightly  affected.  At 
the  very  highest  angles  (near  90  degrees),  the  reflection  coefficients  for  vertical  and 
horizontal  polarization  become  nearly  equal  (exactly  equal  at  90  degrees  since  at  this 
angle  —  vertical  incidence  —  there  is  no  physical  distinction  between  the  polarizations). 


RANGE*  NAUTICAL  MILES 

F  •  1000  HUE  .  SNT  HT  •  80  FT  *  VERT  8H  *  10  DEC  • 

WAVE  HT  .  2  FT  «  BM  TILT  ■=  0  DEC  •  FS  RN0  «  SO  NN  ■ 

PQiRR JZRT ION*  HORIZONTAL 

Fig.  11  -  Detection  oontour  plot  for  the  same  parameters 
as  Fig.  10  except  sea  wave  height  increased  to  2  ft 
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F  =  1000  MHZ  «  ANT  NT  -  80  FT  .  VERT  BW  =  10  DEG  ■ 

WAVE  HT  «  0  FT  ■  8M  TILT  ■>  0  DEG  *  FS  RNG  •  50  NN  « 
POLARIZATION,  VERTICAL 

Fig.  12  -  Detection  contour  plot  for  the  same  parameters 
as  Fig.  10  except  vertical  polarization 


The  last  of  this  series  of  plots,  Fig.  13,  is  the  only  one  for  which  the  antenna  beam 
is  tilted.  The  tilt  is  3  degrees.  This  causes  f  (0)  to  be  different  for  the  direct  and  re¬ 
flected  rays,  and  thus  produces  strange  effects  on  the  interference  pattern,  especially  in 
the  sidelobe  region.  The  disappearance  of  the  lobe  structure  at  about  8-degree  elevation 
is  the  result  of  the  coincidence  of  the  first  null  of  the  free-space  antenna  pattern  with  the 
direction  of  the  reflected  ray.  (The  first  pattern  null  for  a  10-degree  beamwidth  and  the 
assumed  (sin  x)/x  pattern  shape  occurs  at  11.3  degrees  from  the  beam  maximum,  which 
is  here  elevated  by  3  degrees.) 

These  plots  do  not  show  the  full  contour  of  the  lower  side  of  the  lowest  lobe.  The 
contour  is  terminated  at  an  elevation  angle  slightly  greater  than  zero,  because,  as  dis¬ 
cussed  in  the  preceding  section,  the  ray-optical  calculation  of  F  is  not  valid  in  the  "in¬ 
termediate  region."  For  very-low-altitude  targets,  when  this  part  of  the  lobe  pattern  is 
important,  a  second  type  of  plot  has  been  developed. 


Signal -Strength-vs-Range  Plot  for  Constant-Altitude  Target 

In  this  type  of  plot,  exemplified  by  Fig.  14,  the  method  described  for  calculating  F 
in  the  intermediate  region  has  been  used,  for  the  part  of  the  plot  below  that  at  which 
interference-region  calculations  are  valid,  Kerr’s  "bold  interpolation"  is  shown  as  a 
dashed  extension  of  the  signal-strength  curve.  The  input  parameters  of  the  subroutine, 
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Fig.  13  -  Detection  contour  plot  for  the  same  parameters 
aa  Fig.  10  except  3-degree  beam  tilt  angle 
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Fig.  14  -  Plot  of  signal  strength  relative  to  minimum  detectable  signal ,  as  a 
function  of  range,  for  frequency  3000  MHz,  ani\enna  height  100  ft,  target  height 
200  ft,  and  free-space  detection  range  50  nautical  miles 


named  FRPLOT,  are  (in  addition  to  the  x  and  y  dimensions  of  the  plot  in  inches)  the 
frequency  FMHZ  (megahertz),  the  antenna  height  Hi  arid  the  target  height  H2  (feet),  the 
sea-wave  height  (crest-to-trough,  feet),  the  polar izatioA  (horizontal  or  vertical),  and  the 
calculated  or  assumed  free-space  maximum  range  of  the  radar  RFS  (nautical  miles)  (7). 
The  radar  horizon  is  first  calculated  from  the  equation  g^ven  by  Kerr  (1), 

rh  =  v/2^:  (v/h;  +  yhj)  (4i) 


Rh  =  1.2287  (VF7  +  /R7)  .  (41a; 

in  which  ap  is  the  "effective  earth's  radius,"  here  taken  to  be  4/3  times  the  actual  ra¬ 
dius  a.  For  a  =  3440  nautical  miles,  hs  and  h3  in  feet,  and  Rh  in  nautical  miles,  Eq. 
(41a)  applies.  (In  Eq.  (14),  incidentally,  the  denominator  is  equal  to  Rh.) 

The  x-axis  com  iinate  system  is  then  established  for  a  maximum  range  slightly  be¬ 
yond  the  horizon.  In  the  interference  region,  the  calculation  of  F  is  made  from  Eqs.  (7) 
through  (19).  The  signal-level  computations  are  made,  after  calculating  f  at  range  K, 
from  the  formula 


SdB  r  40  log,0(FR0/R) 


(42) 


where  R0  is  the  free-space  radar  range.  This  means  that  the  zero-decibel  level  corre¬ 
sponds  to  the  minimum  detectable  signal  — that  is,  the  signal  level  in  free  space  at  the 
free-space  maximum  range.  Thus,  portions  of  lobes  lying  above  the  zero-decibel  line 
correspond  to  regions  of  signal  detectability,  and  those  below  to  regions  of  no  detection. 
To  make  a  plot  of  this  type  for  a  one-way  radio  system,  it  is  necessary  only  to  change 
the  factor  40  in  Eq.  (42)  to  20. 

Following  the  recommendation  of  Kerr  (1),  the  interference -region  calculation  of  F 
was  terminated  at  the  point  where  the  path  difference  8  of  Eqs.  (1)  and  (10)  is  equal  to 
A./4.  At  greater  ranges,  the  interpolation  procedure  for  the  intermediate  region  was 
used. 

The  wave  polarization  is  taken  into  account  in  the  interference -region  calculations, 
but  those  for  the  diffraction  region  are,  strictly  speaking,  valid  only  for  horizontal  po¬ 
larization.  However,  for  the  frequency  range  being  considered  the  results  for  horizontal 
and  vertical  polarization  are  not  significantly  different  (Kerr,  Ref.  1,  p.  124).  The  an¬ 
tenna  pattern  is  not  taken  into  account  because  for  low-altitude  targets  (for  which  this 
type  of  plot  is  primarily  intended),  the  significant  radiation  all  comes  from  the  same  part 
of  the  beam  — the  part  directed  toward  the  horizon.  The  free-space  range  R0  (Fortran 
parameter  RFS)  should  be  calculated  for  the  free-space  antenna  gain  in  that  direction, 
which  will  usually  be  the  maximum  gain. 

The  signal-level  plot  is  begun  at  maximum  range  and  is  terminated  when  the  range 
has  decreased  to  one-tenth  of  the  total  horizon  range,  Rh,  of  Eq.  (41). 


PLOTTING  TECHNIQUES 

The  plots  shown  in  Figs.  6  through  14  were  made  by  the  Gerber  Model  875  Automatic 
Drafting  Machine  located  in  the  NRL  Engineering  Services  Division.  This  machine  has  a 
maximum  plotting  area  of  5  by  8  feet  and  a  plotting  accuracy  of  0.005  inch  or  better,  de¬ 
pending  on  the  area  of  the  plot.  Standard  Leroy  India-ink  pens  are  used  in  a  rotatable 
turret  so  that  more  than  one  pen  width  can  be  used  in  a  plot.  (The  computer  plotting  sub¬ 
routines  contain  pen-changing  commands  to  the  plotter.)  The  sample  plots  shown  v.ere 
drawn  to  a  size  of  about  10  by  12  inches,  exclusive  cf  lettering,  and  photographically 
reduced. 

The  lettering  and  numbering  is  done  by  the  plotter.  Subroutines  RHACHT  and 
FRPLOT  contain  the  logic  required  for  generating  the  appropriate  numbering  for  the 
specified  values  of  RMAX  and  HMAX  or  H  l  and  H  2.  The  plotting  and  labeling  of  the 
signal-level  coordinate  system  (Fig.  14)  is  done  by  subroutines  named  GRID  and  AXLABL 
which  are  not  listed  in  Appendix  B  because  they  are  lengthy ’and  not  relevant  to  the  pri¬ 
mary  purpose  of  this  report.  The  paper  tape  that  controls  the  plotter  is  produced  by  a 
Fortran  subroutine  package  developed  at  the  Naval  Research  Laboratory  for  this  pur¬ 
pose;  these  subroutines  are  also  not  listed  in  Appendix  B.  The  plotter  has  the  capability 
of  drawing  a  continuous  or  a  dashed  line,  as  illustrated  in  Fig,  14. 

The  interpolation  required  for  the  dashed-line  extension  of  the  signal -level  curve  in 
Fig.  14  is  done  by  computing  two  points  in  the  interference  region  and  two  points  in  the 
diffraction  region,  then  solving  for  the  cubic  equation  of  a  curve  that  passes  through 
these  four  points.  This  equation  is  then  used  as  the  basis  for  the  interpolation  —  that  is, 
for  plotting  the  dashed  section  of  the  curve.  This  computation  is  done  in  a  Fortran  sub¬ 
routine  named  CURVE,  of  which  a  listing  is  given  in  Appendix  B.  The  subroutine  was 
originally  written  in  a  more  general  form  for  plotting  a  smooth  curve  that  passes  through 
any  given  set  of  points;  that  subroutine  was  named  FCURVE,  an  acronym  for  "French 
curve."  The  CURVE  subroutine  is  a  simplification  of  FCURVE  for  the  special  case  of 
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just  four  points.  A  subroutine  named  MATALG  (matrix  algebra)  is  called  by  CURVE. 
Subroutine  MATALG  solves  any  set  of  simultaneous  linear  equations,  as  required  in  this 
case  for  obtaining  the  coefficients  of  the  desired  cubic  equation.  This  subroutine  was 
obtained  from  the  NRL  Research  Computation  Center  program  library;  it  has  the  Co-op 
Identification  F2  CODA  MATALG. 

A  special  plotting  procedure  was  devised  to  handle  the  discontinuities  of  the  curves 
where  they  go  above  or  below  the  maximum  or  minimum  coordinates  of  the  graphs,  as  in 
Fig.  14.  As  each  point  to  be  plotted  is  computed,  it  is  tested  to  determine  whether  it  lies 
outside  the  plotting  area;  if  it  does,  and  if  the  preceding  point  does  not,  then  the  inter¬ 
section  of  the  straight  line  between  these  points  and  the  bdtmdary  of  the  graph  is  com¬ 
puted,  and  the  plot  is  made  to  that  intersection  point.  The  pen  is  then  lifted,  and  the 
plotting  is  resumed  at  the  next  intersection  of  the  computed  curve  with  the  boundary  of 
the  graph.  The  intersections  are  computed  in  a  Fortran  subroutine  called  INTRSECT, 
which  is  not  listed  in  Appendix  B. 

Subroutines  RHACHT,  LOBES,  and  FRPLOT  contain  calls  to  a  subroutine  named 
MINITAPE,  which  is  also  not  listed  in  Appendix  B.  Its  purpose  is  to  minimize  the  amount 
of  paper  tape  used  to  control  the  Gerber  plotter.  It  does  this  by  testing  each  successive 
computed  point  to  determine  whether  the  curve  being  plotted  has  deviated  significantly 
from  a  straight  line  since  the  last  point  plotted;  if  it  has,  the  plotter  pen  is  moved  to  the 
last  preceding  computed  point,  but  if  not,  the  pen  is  not  moved  (no  paper  tape  punch  is 
made).  This  technique  is  possible  because  the  Gerber  plotter  contains  its  own  interpola¬ 
tion  logic  to  move  the  pen  in  a  straight  line  between  successive  points  designated  by  the 
paper-tape  punches. 

The  computer  time  required  for  these  plots,  using  the  CDC  3800,  is  of  the  order  of  a 
minute  or  two  per  plot,  depending  on  the  number  of  lobes  in  the  plot.  The  plotting  time 
on  the  Gerber  machine,  for  plots  of  the  kind  shown  in  Figs.  0  through  13,  is  about  20 
minutes  per  plot.  The  length  of  the  punched  paper  tape  is  generally  a  few  hundred  feet. 

A  rough  estimate  of  the  computer  and  plotter  charges  for  each  plot  is  $25. 
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Appendix  A 


APPROXIMATE  EQUATIONS  FOR  SPHERICAL-EARTH 
INTERFERENCE  CALCULATION 


The  approximate  formulas  for  the  path  difference  £  and  the  grazing  angle  for 
spherical-earth  reflection  are  derived  as  follows.  The  situation  considered  is  shown  in 
Fig.  Al.  As  indicated,  the  target  is  assumed  to  be  so  distant  that  the  direct  ray  and  the 
reflected  ray  are  virtually  parallel.  The  path  difference  b  is  then  obviously  given  by 

b  -  dj  -  d,  .  (Al) 

where  d2  is  the  distance  from  the  antenna  to  the  reflection  point  and  d,  is  the  distance 
that  the  direct  ray  travels  to  reach  the  point  D.  This  point  and  the  reflection  point  B  axe 
equidistant  from  the  target.  The  dashed  line  BD  is  perpendicular  to  AD.  The  problem 
is  to  evaluate  the  path  difference  s  and  the  angle  0  as  functions  of  the  target  elevation 
angle  0d  and  the  antenna  height  h. 

It  is  readily  deduced  that  0  -  0d  +  y,  where  >  is  the  angle  at  the  earth's  center 
subtended  by  the  antenna  and  the  reflection  point.  This  follows  from  the  fact  that  the 
reflected  ray  is  parallel  to  the  direct  ray;  its  elevation  angle  with  respect  t q  the  horizon¬ 
tal  at  the  antenna  is  therefore  9d.  The  horizontal  at  the  reflection  point  is  tilted  an 
amount  equal  to  y  with  respect  to  the  horizontal  at  the  antenna.  I 


Fig.  Al  -  Befleotlon  diagram  for  spherioal 
earth  and  a  distant  target 


25 


26 


L.  V.  BLAKE 


The  first  problem  is  to  evaluate  y,  which  can  be  done  by  considering  the  triangle 
ABC.  The  angle  p  is  given  by 


p,  =  J  +  y  +  eA  ,  (A2) 

since  CB,  the  earth's  radius,  is  perpendicular  to  the  horizontal  at  point  B,  and  «  =  y  +  od. 
Consequently 


a 


-  2y  -  eA 


(AS) 


Applying  the  law  of  sines  gives 


sin 


(f  ■ ar  - s") 


.  /  w 

smlT4 


y  +  0* 


(A4) 


The  trigonometric  functions  can  be  rewritten  as  follows,  using  the  standard  formula 
for  the  sine  of  the  sum  of  two  angles:  L 


sin(r" 2y  “  0d),=  cos('T~ 2y) 8in  (‘~6^ +  sin(f  - 2>) cos  6* 

-  “Sin  (27)  sin  (0d)  +  cos  (2y)  cos  (0d)  ,  (A6) 

8in^+  y  +  =  cos ^  +  yj  sin  0d  +  sin^-  +  yj  cos  6>d 

=  -sin  (y)  sin  0d  +  cos  (y)  cos  (0d)  .  (A6) 

In  order  to  be  able  to  solve  for  y,  it  is  necessary  to  make  the  following  approxima¬ 
tions  which  are  valid  if  y  «  1: 


sin  (y)  »  y  , 


sin  (2y)  »  2y  , 
cos  (y)  *  1  -  ya/2  , 
cos  (2y)  »  l  -  2y*  . 


(A7) 


Making  these  substitutions  in  Eq,  (A4)  results  in  a  quadratic  equation  for  y,  which  has 
the  solution 


y  =  >/[tan  (Od)/3)2  *  2h/3a,  -  tan  (0d)/3 


(A8) 


It  is  evident  that  when  (tan  ed/3)a  »  2h/3o„,  y  is  a  very  small  angle.  To  evaluate  it 
accurately  for  this  case,  it  is  preferable  to  use  an  approximation  which  avoids  the  diffi¬ 
culty  of  evaluating  a  small  difference  between  two  relatively  large  numbers;  the  approxi¬ 
mation  is  based  on  the  relation: 
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y/TTT  *  1  +  ~  -  —  ,  «  «  1  ,  (A9) 

2  8 

which  is  obtained  by  expanding  the  left-hand  side  In  a  Maclaurln  series  and  retaining  only 
the  first  three  terms.  Thus 


y  -  (tan  8d/3)  y/i  +  6h/(ae  tmi  8d)  -  1  (A8a) 


*e  tan  8d  2a2  tan*  8d  10 

This  expression  becomes  increasingly  accurate  as  ed  increases,  while  the  "exact" 
expression  is  subject  to  increasing  numerical  error  as  8d  increases,  if  a  fixed  number 
of  significant  figures  is  carried  in  computation.  For  h  *  100  ft,  if  as  many  as  ten  signif¬ 
icant  figures  are  carried,  the  exact  expression  is  better  than  the  approximation  up  to 
ed  a  10  degrees,  while  the  reverse  is  true  above  ed  =  20  degrees.  , 

7 

■  When  y  has  been  found,  da  can  be  calculated  by  applying  the  law  of  cosines  to  tri¬ 
angle  ABC,  again  using  the  approximation  cos  y  =  l  -  y2/ 2  The  result  is 


and  therefore 

8  =  d,[l  -  cos  <2*d+  2y)J 
=  2d j  sin*  ( 8d  +  y)  , 

Thus  finally 

8  =  2/h1  ♦  a0(a,  +  h)  ya  sin*  (0d  +  y)  . 


(All) 

(A12) 


(A13) 


(A14) 


It  can  be  shown  that  this  equation  becomes  asymptotically  equal  to  the  flat -earth 
expression,  Eq.  (23),  as  y  -*o,  by  the  following  reasoning.  The  ground  range  G  from  the 
antenna  to  the  reflection  point,  Fig.  Al,  is  given  by 

( 

G  s  yae  .  (A15) 


Therefore,  assuming  h  «  ae> 

h*  +  #e(ac  +  h)  y*  »  h*  +  0*  .  '  (Aid) 


If  the  earth  is  nearly  flat,  then 


h*  +  o*  *  d,* 


(A17) 
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Appendix  B 

LISTINGS  OF  FORTRAN  SUBROUTINES 


In  this  appendix,  the  principal  Fortran  subroutines  and  function  subprograms  used 
for  the  two  lobe-plotting  procedures  are  listed.  Some  auxiliary  subroutines  are  not 
listed  because  they  are  not  relevant  to  the  actual  lobe  plotting. 

The  subroutines  and  subprograms  listed  are: 

Subroutine  RKACHT  (1) 

f  Subroutine  ATICK  (1) 

Function  Fl  (1) 

Subroutine  SIMC0N  (1) 

Subroutine  L0BES  (1) 

Subroutine  SEAREF  (1,2) 

Subroutine  FRPL0T  (2) 

Subroutine  INVERT  (2) 

Function  ESS  (2) 

Subroutine  DIFFRACT  (2) 

Subroutine  UFCN  (2) 

Subroutine  CURVE  (2) 

Subroutine  MATALG  (2) 

/ 

The  parenthetic  numerals  indicate  whether  the  subroutine  is  used  in  the  range-height 
angle  coverage-diagram  plots  (1)  or  in  the  constant-target-height  signal-level  plots  (2 y. 
Subroutine  SEAREF  is  used  in  both.  The  purposes  and  logic  of  these  subprograms  are 
described  by  comment  cards.  / 

/ 

The  names  and  brief  descriptions  of  the  subroutines  not  listed  are:  / 

Subroutines  PL0T,  NUMBER,  SYMB0L  (1,2).  These  are  in  the  CDC  3800/on-line 
library  plotting  subroutine  package  for  the  Calcomp  plotter.  Additional  subroutines  with 
entry  points  named  PENCHG,  DASH0N,  and  DASH0FF  are  used  in  Gerber  plotting. 

Subroutine  DEGREE  (1)  plots  a  degree  symbol  for  labeling  the  angi^  scale  of  the 
range-height-angie  chart.  / 

Subroutine.  CENTER  (1,2)  makes  a  computation  required  for  centering  numbers  and 
captions  in  coordinate  labeling.  / 

Subroutine  INTRSECT  (1,2)  computes  the  intersection  point  of  two  straight  lines  when 
the  x-y  coordinates  of  any  pair  of  points  is  given  for  each  line. 


;>W1&TW"'<  '*» '^■"i vr™Tri>ipf*vp.*e**“  '■•"* 
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Subroutine  MINITAPE  (1,2)  checks  each  computed  point  of  a  plot,  and  makes  a  call 
to  PL0T  if  and  only  if  the  point  will  cause  the  line  being  plotted  to  deviate  significantly 
from  a  straight  lino.  This  procedure  minimizes  the  amount  of  paper  tape  required  for 
the  Gerber  plotter  / 

r 

Subroutines  GRID  and  AXLABL  (2); plot  and  label  the  axes  of  the  coordinate  grid  used 
in  the  constant-target-height  signal -l^fel  plots. 


Subroutine  ARR0W  (2)  plots  the  arhow  that  designates  the  radar  horizon  line  in  the 
signal-level  plots.  :  1 

Listings  of  these  subroutines' can  be)  furnished  to  Government  agencies  and  their 
contractors  upon  request.  Requests  to  the  author  may  be  addressed  as  follows:  Code 
5370,  Naval  Research  Laboratory,  Washington,  D.C.,  20390. 


I 

i 
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subroutine  rhachtixmax.ymax.rmax.hmax.anthitE) 
external  FI 

DIMENSION  SN( 181 ) «  RNGK181).  TNI (181).  JA(6) 

DIMENSION  XZ( 180) «YZ( 180) 

DIMENSION  IRFAC  <5>t  IHFAC  <8) 

COMMON  /A/  REF.  GRAD.  RAD*  CONST.  U<182).  N 
COMMON  /MTPE/  X22. Y22 . XI 1 . Y I  1 . XAA . YAA ♦ ERROfe 
DIMENSION  I ANG < 9 ) . NANG ( 9 ) 

COMMON  /&/  XXI 131).  YY (181).  CT1I181).  SNl(lol).  DEL 
COMMON  /XCYC/  XCOR.YCOR 
DATA  (ERROR  =  .001) 

DATA  < ( !ANG< IN) . IN*1 .  9>*l.ll.  31.  5 1 . 1 0 1 « 1 1 1 . 1 21 « 1S1 . 18 1 ) 

DAT  A ( <  NANG  < I G ) . I G= 1 ♦  9)*0.t.  3.  ,,5.10.20.30.60.90) 

DATA  <REF=. 000313) ♦ ( GRAD® . 00004384822 ) 

DATA! < I HFAC (  I  ) . 1  =  1 . 8 >  =50 . 1 00 ♦ 500 . 1000.5000. 10000.50000. lOOOOO) 
DATA  <(  IRFAC  (I).  I  =  1.  5)  =  5.10.50  .100.500) 

DATA  ((JA(IA).  IA  =  1.  6)  =  11.  31.  51.  101.  121.  151)  - — * 

C  FOR  PLOTTING  RANGE -HE 1GHT-ANGLE  CHART  ON  GERBER  PLOTTER 

C  START  WITH  HEAVIER  OF  TWO  PENS.  IN  TURRET  POSITION  11 

CALL  PENCHGI 1 1 ) 

DEL*XMAX# .0 1 

IF  ( YMAX.LT. XMAX )  DEL=YMAX#.01 

BA  *  6076.1155*RMAX*YMAX/<HMAX*XMAX> 

BA2  =  BA#BA 

PREC  ■  .iOOOOl 

CONST  *  .3048/' 1852.0 

RAD  «  20898950.13  +  ANTHITE 

AB  »  1.0  +  REF 

AB2  =  A8#AB 

CD  «  2.0  *  REF  ♦  REF#REF 
I  I  «  0 

DO  29  J I  =  I  .  2 
IF  <JI  «EO.  1)  49.  50 

49  M  «  1 

MM  s  tOO 
GO  TO  51 

50  M  =  11 
MM  «  91 

51  DO  29  IK  *  M.  MM 

IL  ■  < IK-I)*10*#{ JI  -  1  ) 

I  I  =  I  I  +  1 

ELEV1  «  IL 

ELEV  «  El.EV  1/10.0 

RDN  »  ELEV./57. 29577957 

SN( II)  «  SIN(RDN) 

S  •  SN<  I  I )*#2 
U < IJ)  *  AB2#S-CD 
IF  (It  .LT.  181)  360.  361 
360  TN  •  TANF(RDN) 

RDN I  •  ATANF ( BA#TN ) 

TNI ( I  1 )  ■  TANFtRONl ) 

CT 1 (II)  ■  COS ( RDN 1 ) 

SN1 (II)  •  S I N  <  RDN 1 ) 

GO  TO  29  t 

36)  CT1 Cll)  ■  0.0 
SNl ( II)  ■  1.0 
29  CONTINUE 
B=BA#XMAX 
DO  290  IX=  1.1 80 


;  :\)  - 

. . / _ 
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290 


46 


47 


IF 
IF 
I J 
GO 
I J 
OO 
H2 
DO 


XZ( IX)  =  XMAX/(SORT<  1.0  +  ( ( TNI ( IX) >*f2)/BA2> > 
YZ(IX)  *  B*SQRT ( 1 • 0—  C  XZ  < IX ) /XMAX)  **2  >y 
COMPUTE  POINTS  ON  CONSTANT-HEIGHT  CURVES 
HI  =  0.0 
IJM  =  100 

UMAX  ■  |NTF( ALOGIO(HMAX)  -  .477122] 

HINT  *  10.*#( JMAX-1 )  -  l.E-9 
I  JO  =  (HMAX/( 10.**( JMAX-1 > )+  .000(60001) 

DO  30  J  =  It  JMAX 

<J  .EQ.  JMAX)  IJM  a  I  JO 
(J  .EQ.  1)  2.  3 
=  10 
TO  4 
a  20 

30  I  *  IJ.IJMtlO 
*  I  #  104*  ( J  -  1 > 

304  K  3  1,  lei 
N=182-K 

IF  <H2  .EQ.  10  .)  RNGl(N)  /*  0.0 
IF  (HZ  .EQ.  10  .  .AND.  N  AEQ.  1 )  6t  7 
GAM  =  REF*GRAD/AB  j 

Rt  *  1.0/RAD  / 

QO  *  2.0*(RI  -  GAM)  / 

RNG2  *  <C0NST*AB/QQ)*2/.0*SQRT<QQ*H2) 

GO  TO  8  / 

CALL  SIMCON  (H1.H2.PREC  . 1 5 .RI NC t NOl t R.F1 > 
RNG2  =  RNG1(N)  +  RINC 
RNG1(N)  =  RNG2 
IF  <H2  .LT.  HINT)  (SO  TO  304 
A  =  RNG2#XMAX/RMAX 
B  =  BA* A  / 

AH  «  A*A  / 

B2  a  B*B  j 

IF  < N  .EQ.  181V  46.  47 
CALL  PLOT  (O.^e.3) 

XX( 181 )  =  O. 

YY< 181 )  ■  8  / 

XsO. 

Y  =  B 
X22=0« 

Y  2?. 

GO  TO  304 
XLAST 
YLAST  Y 

X  «  A/(SQRT< 1 .0+( ( TNI <N  ))**2)/BA^»» 

Y  •  B^SQRT  (1.0  -  X*X/A2> 

IF  (HZ  .EQ.  HMAX)  87.  68 


2)  780.  781 


XZCN)  +  .0001)  782.  783 
INTRSECT  (XLAST. YLAST. X.Y.  XZ ( N ) • YZ ( N ) * XZ ( N4 1) . VZ < N+ 1 ) «X0  .  YO > 
MINI  TAPE (XO.YO) 

PLOT  ( XO.YO. 2 ) 


/ 


n  n 
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IF  (H2.EQ.HMAX)  785.305 
783  IF  <K  .EO.  181)  787*  788 

787  CALL  MINlTAPE(X.Y) 

CALL  PLOT (X. Y.2) 

XCOR  *  X 

VC OR  «  Y 
GO  TO  304 

788  CALL  M  I N I  T APE  <  X  «  Y ) 

304  CONTINUE 

305  HI  »  H2 
30  CONTINUE 

GO  TO  789 

785  CALL  PLOT ( XZ (l>»0»*3> 

DO  786  NX  =  1  .  N 
XX(NK) =XZ(NK> 

YY<NK)=YZ(NK) 

786  CALL  PLOT ( XX  <  NX ) < YY ( NX )  « 2) 

CALL  PLOT  <  XO  «  YO . 2 1 

XCOR  a  XO 
YCOR  ■  YO 

C  DRAW  X-Y  AXES 

789  X  =  0. 

Y  *  YMAX 
CALL  PL0T(X«Y.3) 

Y  »  0.0 

CALL  PLOT <  X. Y.2  > 

XfXMAX 

CALL  PLOT (  X. Y  .2  ) 

DRAW  CONSTANT-RANGE  ELLIPSES 
CHANGE  TO  FINER  PEN 
CALL  PENCHGI 10) 

KAM=RMAX-1 . 

INTR=|0 

IF  (RMAX.LT .  100.  >  INTR>=5 

IF ( RMAX.GT .  300.)  I NTR»25 
DO  31  KA= INTR.KAM, INTR 
RG  o  XA 

A  *  RG*XMAX/RMAX 
B  s  BA*A 
A2  -  A#A 

DO  38  XC  *  1*  181 

XR=XC 

IF  (KC  .EO,  181)  91*  92 
91  CALL  M I N I T APE ( O.«0 ) 

CALL  PLOT ( 0. <  B«  2 ) 

GO  TO  31 

92  X  m  A/(S0RT(  1  ,0+(  (TNI  (KC>  )##2>/'BA2>  1 

Y  »  8*5GRT  (1.0  -  X*X/A2) 

IF  <  Y  .GT.  (YY(KC)  +  .0001))  73*  72 

73  IF  (KR  .EO.  N+l)  730.  731 

730  XlaXCOR 
Y l a YCOR 
GO  TO  732 

731  Xl»XX(KR-l) 

Y  1  »  Y Y ( KR  -  1  ) 

732  X2  »XX <  XR 1 
Y2  »  YY(KR) 

XBsX 
V8  «  Y 


u 
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l* 

, vr4**  **'  1  * 


I 


i 


CALL  I N TR3KC T  (  X  l  «  Y 1  *X2»Y2»XA,YAt>XO»Y8*XO*YO) 
CALL  PLOT < XC • VO  *2  > 

00  VO  a l 

7a  XA  »  X 
YA  a  Y 

ip  <kc  «kq«  n  n»*  ©6 
as  call  plot <  x, v « 3 i 
xaa-x 

Y8R-Y 

00  to  aa 

as  IF(KC.EQ,2)  860,861 
060  XI  1<*X 
Y||«Y 
XAA«X 
YAA-Y 

00  TO  38 

861  CALL  M IN  I T APE  <  X,  Y) 

38  CONTINUE 
31  CONTINUE 

c  ORAW  ANGLE  TICK  MARKS 

700  CALL  ATICK  11*120,1,4.5) 

CALL  ATICK  <121,183,5.1.2) 

C  DRAW  RAY  LINES 

DO  34  KK  «  1 ,  6 
NP  ■  J  A  <  KP'  1 
X  •  0,0 

Y  «  0.0 

CALL  PLOT  <  X. Y,3 ) 

X  «  XX INF) 

Y  «*  YY  <  NF 1 

CALL  PLOT<X.Y,2) 

34  CONTINUE 

C  MAKE  RANGE  TICK  MARKS  ON  Y  >  0  AXIS 

Kt  ■  9 

DO  364  KY  a  1 ♦  800 
KT  »  KT  +  l 

IF  <KT  »EQ.  10)  368.  369 

368  FAC  »  2.0 
KT  n  0 

GO  TO  370 

369  FAC  =»  1.0 

370  R2  a  ICY  -  1 

X  =  R2*XMAX/RMAX 

IF  ■  <  X  .GT.  <  XMAX  +  .0001)'  GO  TO  365 

Y  r  0.0 

CALL  PLOT  <X. Y.3) 

Y  *  -  FAC  *  DEL 
CALL  PLOT  <  X*  Y ♦ 2 ) 

364  CONTINUE 

C  MAKE  HEIGHT  TICK  MARKS  ON  X  »  P  AXIS 

365  KS  a  9 

K.JM  e  HMAX  /  HINT  +  1.001 
DO  37  KJ  »  I .  KJM 
KS  «  KS  +  l 

IF  ( KS  *EQ.  10)  375,  376 

375  FAC  a  2,0 

KS  =  0 
GO  TO  377 

376  FAC  «  1.0 


o  o 
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377  H  *  HINT  »  t KJ  -1 ) 

V . «  H4YMAX/HMAX 
X  ■  0.0  1 

CAUL  PLOT ()i*  Y,3)  . 

X  - - FAC  #  DEL 

CALL  PLOT ( X*  Y  »2 ) 

37  CONTINUE  * 

LETTER  AND  NUMBER  RANGE «  HEIGHT*  AND  ANGLE  SCALES 
CHANGE  TO  HEAVIER  PEN 

call  penchgi in 

IF  (XMAX-YMAX)  460*460*461 
460  SFAC«XMAX  *.125 
30  TO  462  ’ 

46(1  SFAC-YMAX4.  J2S 
462  M«.I75*SFAc‘ 

C  SELECT  SIZE^  OF  RANGE  AND  HEIGHT  NUMBERING  I  INTERVALS 

OO  100  IR  »'  1  *  5 
NR  »  RMAX  /  IRFAC  (IR) 

IF  (NR  «LE.  10)  «01.  100 

101  I  RUN  IT  «  IRlUc  (IR) 

GO  TO  102  , 

100  CONTINUE 

I RUN I T  « I RF AC ( 5 ) 

102  OO  1 10  IH  a  1  •  6 

NH  s  HMAX  /  IHFAC  (IH) 

IF  (NH  ,L£.  10)  |63«  110 
S  03  I HUN I T  »  IHFAC  ( IH) 

GO  TO  128 
110  CONTINUE 

IHUNIT= IHFAC(B) 

128  X  —  .054SFAC 
Y«-.S4SFAC 

CALL  NUMBER (X«Y*H«0*0.0* 2H II) 

N  *  IRUNIT 

120  IDIGITS  «  ALOG10(FLOATF(N) )  +  1.000001 
CALL  CENTER  (H,  ID)GITS*4*BI AS) 

X  «  (NXRMAX)  »  XMAX  -  BIAS 

IF  (X  +  BIAS  *GT.  XMAX)  GO  TO  801 

CALL  NUMBER  (X*Y*H*N«0.0«  2H 14) 

N  ■  N  t-  IRUNIT 
GO  TO  120 
801  V«— 1 »0#SFAC 

CALL  CENTER  (H«  21*21«B1AS) 

X  «  0,5  «  XMAX  -  BIAS 

CALL  SYMBOL! X.Y.H.21HRANGE*  NAUTICAL  MILES. 0*0*21 ) 

X«~.5*SFAC 

Y«-c0875#SFAC 

CALL  NUMBER ( X«Y *H« 0* 0*0 *2HI I ) 

N  ■  IHUNIT 
X«-t , 1 4SFAC 

121  Y* (NXHMAXI4YMAX  -  .O075#SPAC 
IF  (Y  «GT •  YMAX )  GO  TO  803 
CALL  NUMOER  (X.V.H.N*  0»0*2HI6  ) 

N  ■  N  +  IHUNIT 

GO  70  121 

603  X  «  -1.404SFAC 

CALL  CENTER! H»  12* 12 *01  AS) 

Y  »  0,5  4  YMAX  -  BIAS 

CALL  SYMBOL  (X.V*H» 12HHEIGHT,  FEET, ©0**12) 
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XPR  =  2.*XMAX 
YPR»-2.0*H 
00  804  IL  »  1 ♦  9 
IND  ■  !ANG<  tl_> 

NAN  =  NANGIIL) 

CX=«2 

IF  <IL  »GE«  5)  CX» *  1 25 

XaXXCIND)  +  »4#CT 1 < I  NO ) *SFAC  -  CX*SFAC 
Y=YY< IND)  +  .4*SN1 < INQ)*SFAC  -  .0875»SFAC 
IF  ( IL  *EQ.  9)  X  ■  X  -  H 

IF  UL.EO.  9  «0R«  < Y-YPR)  .GT.  <i»5*H>>  GO  TO  888 
IF  ((XPR-X)  *LT«  <3.0*H>  .OR.  X  .LT.  2<*H>  GO  TO  804 

088  CALL  NUMBER ( X * Y « M « NAN . 0 . 0  «  2H  i  2  1  - - 

1  CALL  DEGREE  <X  +  ,35*SFAC.  Y  +  .175  *  SFAC* .0e«SFAC> 
XPR*X 
YPRsY 

804  CONTINUE 
END 


SUBROUTINE  ATICK  < IA.JA.KA. MB .MCI 

THIS  SUBROUTINE  PLOTS  ANGLE  TICK  MARKS  ON  RANGE-HEIGHT-ANGLE  CHART. 
COMMON  /&/  XX  <  181  )  «  YY  (181).  CT1  <1811.  SNKIBl).  DEL 
MA  b  MB 

DO  1  K  «  IA.JA.KA 
MA  .  MA  +  | 

IF  <  MA  .EQ.  MC)  2.  3 

2  FAC  *  2.0 
MA  ■  0 

GO  TO  4 

3  FAC  *  1 .0 

4  X  ■  XX<K> 

Y  «  YY ( K  > 

CALL  PLOT (X. Y.3 ) 

X  *  X  +  DEL*FAC*CT1 <K) 

Y  «  Y  +  DEL*FAC#SN1 <K) 

CALL  PLOT ( X. Y . 2 ) 

1  CONTINUE 
END 


FUNCTION  FI <  X J 

INTEGRANO  FOR  Si Me ON  SUBROUTINE 

COMMON  /A/  REF.  GRAD.  RAO.  CONST*  U<182>.  N 

BB»REF*EXP( -GRAO#X 1 

CC*X/RAO 

V ■ 2  »  0  #B B+£3*BB 

W*2.0#CC+CC#CC 

FX  «  SORT  (U(N)  V  +  W  ♦  V*W) 

FI  *  CONST  *  <1.0  -*V>*<  1  .0  ♦  CC 1 /FX 
END 


oo 
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SUBROUTINE  SIMCONCXl « XENO tTEST .LIM. AREA.NOt *R«F> 

D1  UCSD  SIMCON  .  ' 

REVISED  OCTOBER  196?  LATER  REVISION  BY  URM  'FEB  1967 

NOIeO 

C  THIS  STATEMENT*  WHICH  READ  #R1*10*»  REMOVED  BECAUSE  NOT  NEEDED* 

ODD-0,0 
I  NT**  1 
V«1  ,0 
EVEN*0.0 
AREA1 *0*0 

19  ENDS*F (XI )+F ( XENO  > 

2  HMXEN0-X1  )/V 
ODb*EVEN+OOD 
X=Xl+H/2. 

EVEN*0.0 

00  3  I  *  1  *  I NT 

21  eV,EN*EVEN+F<X)  '  > 

X=»X+H  , 

3  CONTINUE 

31  AREA*(ENDS+4.0*EVEN+2.0*ODD>«H/6*0 

NU1=*NU1  +  1  A 

34  R=ABSF( (AREA  1 -AREA) /AREA) 

IF ( NO  I  — L IM>  34 1 « 3S«  33  '  x  1  '  '  • 

341  IF  (R-TEST  )3S«3S.«4  ’ 

35  RETURN  '  \  > 

4  AREAJbARFA  \  ’  '  l  \V 

46  I  NT**  24 1  NT  \ 

V»2«0#V.  ,  '  .  •  .. 

GO  TO  2  v  .  >.  ■  " 

END  K  ‘  / 
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SIMC0N22 
SIMC0N23 
S1MC0N24 
SIMC0N25 
SIMCON?* 
SIMCON27 
SIMCON20 

'  SIMCON29 
SIMC0N30 


■M  :  '  ' 


38 


L.  V.  BLAKE 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


SUBROUTINE  LOSES  ( XM AX  . YM Ax  <  RM  A X • HM AX «  FMHZ  «  AHFT  «  B  WD < WHF  T • T I Lt «  THM  « 

■  I  I  POL  *  RF  S ) 

TH!£  SUBROUTINE  PLOTS  SEA-KEFLECT I  ON  INTERFERENCE  HAT TERN  FOR  RADAR  AT 
FMHZ*  ANTENNA  HEIGHT  AHFT  ABOVE  SEA  WATER.  PLOT  IS  MADE  ON 
RANGE -HEIGHT-ANGLE  CHART  OBTAINED  BY  CALLING  SUBROUTINE  RHACHT  WITH  MAXIMUM 

T rHTALl^«f ;  VMAX  INCHES.  MAXIMUM  RANGE  OF  CHART  RMAX  N.  MILES. 
MAXIMUM  HEIGHT  HMAX  FEET.  BWD  IS  ANTENNA  VERTICAL  BEAMWJDTH.  DEGREES. 

WHFT  IS  ASSUMED  CREST-TO-TROUGH  WAVE  HEIGHT.  FEET.  TILT  IS  ELEVATION  ANGLF  op 
ANTENNA  MAXIMUM.  DEGREES.  THM  IS  MAXIMUM  ELEVATION  ANGLE  TO  WHICH  LOBE 

PLOT  IS  DESIRED.  IPOL-1  FOR  VERTICAL  POLARIZATION.  «2  FOR  HORIZONTAL. 

RFS  IS  ASSUMED  FREE-SPACE  MAXIMUM  RANGE  OF  RADAR  FOR  SPECIFIED  TARGET. 

NAUTICAL  MILES  (MUST  NOT  BE  GREATER  THAN  HALF  OF  RMAX). 

SUBROUTINE  WRITTEN  BY  L.  V.  BLAKE.  CODE  S370.  NRL. 1968.  FOR  HORIZONTAL 
POLARIZATION.  MODIFIED  SEPT.  1969  TO  ALLOW  VERTICAL  POLARIZATION. 

,'COMMON  ,'MTp£/  X2»  Y2«  XU.  YU.  XA.  YA.  ERROR 

DATA  ( P  |.*3«  14 1 592654 )  .  < HALFPI « |  «S7p 796327 )  .  (P 1 2«6.283 1 B5307 )  . 

1  (RDN».G1T4S329252) ♦ <COHV« 1 .645788333E-4)  .(AE  »  2.786S266B4E7) 

COMPUTE  ELLIPtICITY.  L«  OF  CONSTANT-RANGE  CONTOURS.  AND  STANDARD 
-DEVIATION,  OF  SINUSOIDAL  WAVE  F  CREST-TO-TROUGH  HEIGHT.  WM'T. 

\  '■  E-  *  VMfX#RMAtfX<XMAX*HNAX#CONV  .  ' 

;  H  -  >,WHrTf  ..&5&S534  .  '  'V  •.  '■  4  :  •  , 

COMPUTE  RADIAN  VALUES  OF  ANGLES  GI\VEN  IN  DEGREES.  , 

TIL'TR  *  TILT  -ft  RON  i  \ 

\  .  THMAX  *  THM1  A'j^bN  .  \  .  >  • 

i  BWR  =  BWO  *'  RON 

CONSt.NT^eOEDN,N^XP«ESSI«-  FOB  CS1N  Xl/X  MW. 

COMPUTE  WAVELENGTH,  W.  *  ' 

' W  *  983.573/FMHZ 
W4  ».  0.25  *  W 

WLIM  «  0.0 I  *  W  I 

FAC  PI  2/f  \  ’ 

COMPUTE  ANGLE  INCREMENT  FOR  PLOTTING  LOBES.' 

DELI  »  A,92/(FMHZ*AHFT)  . 

AH2  «  AHFT  #*  2  ■  j  '  . _ 

HAE  »  2.  #  AHFTV< 3,*Ae> 

AEH  «  AE  *  <  AE  +  AHFT )  i 

PARAM  ■  SORT  (AE/’IEv  *AHFT>  :  ,  \ 

TUPIW  k  Plg/W 

INITIALIZE  TARGET  ELEVATION  ANGLE.  THET2. 

THET2  -  -DELI  \ 

N  *  U  \ 

INDEX  «  0  , 

1  THET2  «  TMET2  +.UELI  x 

IF  <TMET2  «GE*  HALFPI)  GO  TO  40  ' 

T2  *  TANF ( THCT2 ) . 

52  *  StNlTiHETgl 

53  *  S2  • 

PS!  *  THET2  , 

COMPUTE  PATH  DIFFERENCE.  PO.  USING  FLAT-EARTH  FORMULA. 

PD  •  2,  #  AHFT.  *  S2. 

IF  UNDEX  oEO*  t>  GO  TO  77 
T23  «  T2/  3. 

eSthIE  °BAZ,NG  ANGLE*  pS«.  AND  PATH  DIFFERENCE.  PD  I .  FOR  SPHERICAL 
GAM  ■  SORT  (T23*»2+HAE)-T23 

PS  I  «  THET2  +  GAM  „  -- 

S3  «  SINIPSI) 

ZETA  «  PARAM  #  T2 


i 
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DI  3  0.57735  *  SORT!  1.0+2. 0*ZE TA/SURT < ZET A**2+3.  >  ) 

PD  1  =  SQRT< AH2+AFH#GAM4*2) *2.*S3**2 
IF  ( PD  1  «LT.  W4)  GO  TO  1 
N  =  N+  1 

IF  ( ABS ( PD-PD I )  *LT « WL IM ) 78 *  79 

C  WHEN  FLAT-EARTH*  PD.  ALMOST  EQUALS  PDl*  THEREAFTER  L)0  FLAT-tARTH 
c  APPROXIMATION  FOR  PD*  PS  I  *  AND  D  (INDEX  =  1  CAUSES  OMISSION  OF 
C  SPHERICAL-FARTH  CALCULATIONS). 

73  IF  <01  .GE.  .999)  780.  79 
780  Di  =  1. 

INDEX  *  1 
79  PD  s  PD  I 

C  SUBROUTINE  SEARCH  COMPUTES  SMOOTH-SEA  COMPLEX  COEFFICIENT  OF 
C  REFLECTION,  RHO*  PHI,  FOR  GRAZING  ANGLE  PGl  AND  FREQUENCY  FMHZ • 

77  CALL  SFAREF  ( FMHZ *  PS  I  *  I  POL  * RHO *  PH  I  ) 

C  COMPUTE  ROUGHNESS  FACTOR.  REFER  TO  NRL  REPORT  6930*  EG*  S6» 

RUF  =  EXP(-a.*<PI*H*S3/W)**2) 

C  ANG  IS  ANGLL  THAT  DIRECT  .RAY  MAKES  WITH  DIRECTION  OF  UEAi*l  MAXIMUM. 
ANG  =  THET2  -  TILTR 
IF  (ANG  .EQ.  O.)  22*  23 

C  COMPUTE  PATTERN  FACTOR*  PAT.  FOR  DIRECT  W^Y. 

2?  PAT  =  1  . 

GO  TO  24 

23  UU  =  CONST  *  SIN  (ANG) 

I  NT  a  UU/PI2 

DIFF  =  UU  -  I  NT  *  PI? 

PAT  =  SIN(DIFF)  /  UU 

C  ANGR  IS  ANGLL  MAUL  BY  REFLECTED  RAY  WITH  DIRECTION  OF  bfcA|v|  MAXIMUM* 

■  ?4-  ANGR  n  THET2  +  TILTR  +  2.  *  GAM 

IF  (ANGR  «EO.  0.)  25*  26 

25  PATR  =  I . 

GO  TO  27 

C  COMPUTE  PATTLRN  FACTOR,  PATR.  FOR  REFLECTED  RAY. 

26  IJUR  =  CONST  *  SIN  (ANGR) 

INTR  n  UUR  /  PI  2 

niFFR  s  UUR  -  INTR  *  PI2 
PATR  s  SIN(DIFFR)  /  UUR 
urs — (  AQS(PAT  )  .LT.1  ,E-44  )  28*29 
?A  IF  (PAT.LT.O.)  ?frug81 

280  O  3  -RhO*RUF#PaTR*1  .p4B*r>l 

GO  TO  30  X 

281  D3ni*WHO#RUF*PATR*l .645 

GO  TO  30  X 

C  COMPUTE  tFFcCTIVL  REFLECTION  COEFFICIENT.  REFER  TO  NRL  REPORT  6930* 
C  EOLATION  62. 

n  s  DI#  RHO  #  PATR/PAT  #  RUF 

C'  COMPUTE  PHASE  O IFFERENCE *,  ALPHA .  AND  PATTERN-PROPAGAT l ON  FACTOR. F. 

30  ALPHA  n  FAC.  *  PD  +  PHI 
INTI  3  ALPHA/P  1 2 
DIFF1  .=  ALPHA-  INT1#PI2  \ 

F  b  AbS(PAT)#  SORT! 1.0  +  D*D  2. 0*D*CUS ( DI FF I ) ) 

C  COMPUTE  X  AND  Y  COORDINATES. CORRESPONDING  TO  RADAR  RANGE  (KFS#F)  AND 
C  ELEVATION  ANGLE*  THET2*  MOVE  PEN  TO,  POINT  (X*Y). 

A  «  RFS#F#XMAX/RMAX  ’* 

X  b  A/SORT,<  I  ,+T2#T2  ) 

Y  b  F  *  SORT  <A#A-X#X) 

IF  (N.FO.t)  2,3 
2  CALL  PLOT  <X*Y,3) 

X2bX 
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\ 

\  Y2*Y 
YP*Y 
XP»X 
<30  TO  1 

3  IF  <N  tF.Qt  2)  44 «  45 

44  XA«X 
’  YA  =  Y 

YP«Y 
XP*X 

Xll«x 

VI  1=Y 
GO  TO  1 

C  IF  LOBE  GOES  ABOVE  TOP  OF  CHART «  TRUNCATE  IT  AND  DfiAW  A  DASHED  LINE 
C  AT  Y  »  YMAX. 

45  IF  ( Y  .GT.  YMAX)  5.  6  \ 

5  IF  <YP  «LT .  YMAX)  50.  51 

50  CALL  INTRSECT<0.. YMAX. XMAX. YMAX. XP.YP,X.Y»XOf.YO) 

CALL  PLOT(XO. YMAX.2)  v  ... 

Ft  YP  =  Y 

XPeX  ' 

'  GO  TO  1 

6  IF  <YP  .GT.  YMAX)  52.  53 

52  CALL  I NTRSECT  <  p . • YMAX . XMAX . YMAX  «XP. YP. X .Y.XO.YO) 

CALL  PEnCHG  t  4 ) 

CALL  PLOT (XO. YMAX. 2 ) 

CALL  PLOT <X6; YMAX. 3) 

X2  =  X0 
Y2eYMAX 
YA  =  Y 
XA  =  X 
Y1I«YA 
XI 1 »XA 
xp=y 

YPbY 
GO  TO  I 

53  CALL  MtNITAPF(X.Y)  .  ,  — ■ 

C  TERMINATE  PLOTTING  WHEN  ELEVATION  ANGLE  EXCEEDS  THMAX • 

IF  <  THET2  »GE.  THMAX)  GO  TO  40 
XP  a  X 
YP  =  Y 

GO  Vo  i 

41.  CA^-L  PLOT  (  X.  Y.2  ) 

END 
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SUBROUTINE  SEAREF  (FMHZ.  PSl.  IPOL.  RHO.  PHI  I 
C  THIS  SUBROUTINE  COMPUTES  COMPLEX  REFLECTION  COEFFICIENT  OF  SEA  WATER.  AS 
C  FUNCTION  OF  FREQUENCY  FMHZ  IN  MEGAHERTZ.  GRAZING  ANGLE  PS»  IN  RADIANS.  WAVE 
C  POLARIZATION  IPOL  (J  FOR  VERTICAL*  2  FQR  HORIZONTAL).  OUTPUT  IS  MAGNITUDE  RHO 
C  AND  PHASE  ANGLE  PHI  (RADIANS)  OF  COMPLEX  COEFFICIENT.  COMPUTATION  IS  BASED 
C  ON  EOS.  <  l‘)~  AND  12)  OF  RAO  LAB  VOL*  13.  PAGE  396.  SUBROUTINE  WRITTEN 
C  BY  L.  V.  BLAKE  NRL  COOE  S370  SEPT  1969. 

TYPE  COMPLEX  EPSC.  GAM.  SQTERM.  TERM 
DATA  (FLAST  ■  0.) 

SINPSt  •  SIN  (PSD 

CSPSI  a  COS (PS I )##2 

IF  (FMHZ  .EQ'<  FLAST)  GO  TO  200 

C  IF  SUBROUTINE  HAS  BEEN  CALLED  PREVIOUSLY  DURING  PROGRAM.  AND  FREQUENCY  IS  SAME 
C  AS  ON  LAST  PREVIOUS  CALL.  PART  OF  COMPUTATION  NEED  NOT  BE  DONE  SINCE  REQUIRED 
C  VALUES  HAVE  BEEN  STORED. 

FLAST  ■  FMHZ 
C  COMPUTE  WAVELENGTH 

W  ■  299.793  ✓  FMHZ  ' 

IF  (FMHZ  .LE.  1500.)  ISO.,. IS! 

C  SIG  IS  CONDUCTIVITY.  MHO/METER.  EPS1  IS  REAL  PART  OF  DIELECTRIC  CONSTANT.  BOTH 
C  DEPENDENT  ON  FMHZ. 

150  SIG  ■  4.3 
EPS!  >  80. 

GO  TO  IBS 

151  SIG  *  4.3  +  (FMHZ  -  1500. >#  .00148 
,IF  (FMHZ  .L£a  3000.)  1S3.  1S4 

153  EPS I  a  aO.  -  (FMHZ  -  1300. )  *  .00733 
GO  TO  IBS  ‘ 

154  EPS  1  •  69.  -(FMHZ  -  3000.)  #  .002429  ' 

SIG  ■  6.32  +  (FMHZ  -  3000.)  #  *001314 

~1S5  EPSC  a  CMPLX  (EPS1 «-60.*W»SlG> 

200  SQTERMa  CSORT (EPSC -CSPSI ) 

IF  (|POL  «EQ.  1  )  160.  161 

160  TERM  a  EPSC  #  SINPSl 

GAM  a  ( TERM-SQTERM ) ✓ ( TERM+ SQTERM ) 

GO  TO  180 

161  GAM  a  (SINPSl  -  SQTERM)  ✓  (SINPSl  ♦  SQTERM) 

180  RHO  a  CABS  (GAM) 

PHI  a-ATAN2  ( A I MAG ( GAM ) «  REAL  (GAM)) 

ENO 


/ 


/ 


— ‘V 
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SUBROUTINE  FRPLOT  (XMAX. YMAX.Ht .H2.WHFT.FMHZ. IPOL.RFS) 

C  THIS  SUBROUTINE  PLOTS  RADAR  RECEIVED-SIGNAL  LEVEL  IN  DECIBELS  ABOVE 
C  MINIMUM  DETECTABLE*  WITH  SCALE  FROM  -20  TO  +60*  FOR  POINT  TARGET  AT 
C  CONSTANT  HEIGHT*  H2*  ABOVE  SEA.  RADAR  ANTENNA  HEIGHT  IS  HI*  S^A  WAVE 
C  CREST-TO-TROUGH  HEIGHT  WHFT.  FEET.  RADAR  FREQUENCY  FMHZ*  MEGAHERTZ. 

C  I POL  IS  1  FOR  VERTICAL  ANO  2  FOR  HORIZONTAL  POLARIZATION  .  AND  RFS  IS 
C  RADAR  FREE-SPACE  RANGE  IN  NAUTICAL  MILES.  COORDINATE  GR 10  DIMENSIONS 
C  ARE  XMAX.  YMAX.  INCHES.  EARTH  CURVATURE  EFFECTS  COMPUTED  BY  METHOD 
C  DESCRIBED  IN  RAD  LAB  VOL*  13*  PP.  1 I3FF.  (KERR*  PROPAGATION  OF 
C  SHORT  RADIO  WAVES.  MC  GRAW-HILL.  1951).  ROUGH  SEA  EFFECTS  BY  EQ.  56* 
C  NRL  REPORT  6930.  SMOOTH-SEA  REFLECTION  BY  EQS.il )  AND  (2).  P.  396* 

C  KERR.  .VALUES  OF  SKA  COMPLEX  DIELECTRIC  CONSTANT  DISCUSSED  NRL  REPORT 
C  7098.  SUBROUTINE  WRITTEN  BY  L.V.  SLAKE*  CODE  5370*  NAVAL  RESEARCH 
C  LABORATORY.  VERSION  OF  MARCH  1970. 

EXTERNAL  ESS  " 

COMMON/ B/T , TT 

COMMON  /MTP£/  X2» Y2*X I . Y I *XA * YA .ERROR 
COMMON  /MN/  Ml*  N1 
COMMON  /XYL/  XLAST*  YLAST 
DIMENSION  XX ( 4 « 1 ) «  YY<4«1) 

DATA  (ERROR  »  .001) 

OATA  (PI  a  3.1 4 15926S36 ) • ( RDN  a  57.2957795) •( TUPI  «  6.281853072). 

1  (PI2  »  1.5707963268) 

YF(Q)aYMAX#{ .25+.5*AL0G 1 0 ( O ) ) 

INDFX  »  0 
FRAC  ■  .99 
TIC  =  YMAX  #  .01 
AU  *  0.0 - 

C  CONVERT  CREST-TO-TROUGH  WAVEHE I GHT  TO  STANOARD  DEVIATION  OF  SINUSOID. 
WH  a  WHFT  *  .35355339 

C  COMPUTE  W  «  WAVELENGTH.  FEET.  AND  PARAMETER  T*  - 
w  a  9B3.S73/FMHZ  ' > 

T  a  SORT  (H1/H2)  "  • 

IF  (T  .GT.  1.0)  T  a  1.0/T 
TT  a  T  #  T 

C  COMPUTE  RH  a  TOTAL  HORIZON  DISTANCE.  NAUT.  MI. 

RH  *  1.2289  *  (SQRT(Hl)  +  SORT  <H2>> 

W4  *  W  #  .25 

FAC  a  2.0  *  HI  #  H2/6076. 1 155 
HH» ( (H2-H1 )/6076.1 155)#*2 
C  DRAW  COORDINATE  GRID.  LABEL  AXES. 

NCX  -  I NTF  (0.5  #  RH+  l.E-9)  +2 

RMAX  »  2.  *  NCX 

XRaXMAX/RMAX 

XRH  »  XR  *  RH 

RTEST  a  0.1*RMAX 

CALL  GRID  ( XMAX . YMAX. T IC « T IC .NCX.8* 1.1.2.2.1.10.1*10) 

HLETTER  a  ,0173  #  YMAX 

XS  a  XRH  -  12.  #  HLETTER 
YS  a  YMAX  ♦  HLETTgR 

CALL  SYMBOL  (XS.  YS.  HLETTER.  I 3HRAOAR  HORIZON.  Q.O.  13) 

VAR  a  YS  ♦  HLETTER 

YAT  a  YMAX  Q  .  I  #  HLETTER 

TIP  a  HLETTER  *  0.7 

CALL  ARROW  (XRH.  VAR.  XRH.  YAT.  TIP) 

ANTX  a  2, 

CALL  AXLABL  ( XMAX. YMAX. HLETTER. 0. .ANTX. RMAX. O. I »-20. *20. .60. .0. I . 

I  21 HRANGE .  NAUTICAL  MILES. 21 .22HSIGNAL  LEVEL*  0ECIBELS.22I 
CALL  OASMON 
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C  DRAW  FREE-SPACE  SIGNAL  LEVEL.  DASHED. 

DEL  -  XMAX*. 0 1 
X a XR4RFS* ,03 1 6227766 
NN-  <  XMAX— X ) /DEL 
CALL  PLOT  (X.YMAX.3) 

XLAST  *  X 
YLAST  »  Y 

00  1000  12  *  1.  NN 
X  «  X  +  DEL 
R  «  X/XR 

DB1-40.*AL0G10(RFS/R> 

Y - YMAX# ( • 25+ .01 25*0B 1 ) 

IF  (Y.LT.O.)  1001.1002 

1001  CALL  INTRSECT  ( XLAST. YLAST .X .Y . 0. ♦ 0. .XMAX .0. .XO .YO ) 

CALL  PLOT  < XO ♦ YO. 2 ) 

GO  TO  1003 

1002  CALL  PLOT  (X.Y.2) 

XLAST  «  X 

YLAST  -  Y 
10O0  CONTINUE 

C  DRAW  HORIZON  LINE.  DASHED. 

1003  CALL  PLOT  (XRHt  YMAX.  3) 

CALL  PLOT  (XRH.  0,.  2) 

C  START  COMPUTATION  AT  RANGE  20  PERCENT  GREATER  THAN  THAT  AT  WHICH  PATH 
C  DIFFERENCE  IS  QUARTER  WAVELENGTH.  OR  AT  HORIZON.  WHICHEVER  IS  LEAST. 

R  «  FAC/W4  *  1.2 
IF  (R  «GT.  RH)  R  ■  RH  «  .999 
XLAST . *  XMAX 
YLAST  -  l.E-6 
1  IF  (AJ  -  0„999)  4.  4.  5 

5  IF  (D  -  0.9991  4.  4.  6 

C  WHEN  RANGE  HAS  DECREASED  TO  POINT  AT  WHICH  EARTH  CURVATURE  IS 
C  NEGLIGIBLE.  COMPUTE  THEREAFTER  ON  FLAT-EARTH  BASIS. 

6  AJ  *  1.0 
O  ■  I  * 

A«  «  1. 

GO  TO  3 

C  COMPUTE  PARAMETER  S  (FRACTION  OF  TOTAL  HORIZON  RANGE). 

4  S  •  R/RH  1 

SIN  *  S 
SIM  -  1 .2  *  S 

C  FIND  PARAMETER  Si  BY  ITERATION  OF  FUNCTION  ESS.  WHICH  OEFINES  S  AS 
C  FUNCTION  OF  Si, 

CALL  INVERT  (S1N«S1M*4. 1 5. NO I .51 . S .ESS) 

551  -  SI  #  SI 

‘  SOI  -  (1.  -  SSI)  **  2  +  4.  *  SSI  *  TT 
C  COMPUTE  S2*D.J.K  OF  EQS.  ON  P.ltS.  RAO  LAB  V0L.13  (Dl-D.  AJ-J.  AK-K). 
S2  ■  (SORT  (SOI)  -  U  +  SSI)  /(g,  *  SI  *  T) 

552  -  S2#S2 

SQ2  ■  l.+(4.*SSl#S2*T)/(S«(l.-SSl>*(l.+T>)  ( 

01-  1 ./SORT ( SQg ) 

«  ( 1 ,-SSl ) #( 1 «— SS2) 

AK  «  ((l.“SSl)TTT*(|.-SSg))/(l.+TT> 

C  COMPUTE  PATH  DIFFERENCE.  DELTA.  AND  SLANT  RANGE.  R SLANT,  AT  DELTA  « 

C  W4«  START  PLOTTING. 

-  3  DELTA  ■  FAC/R  #  AJ 

RSLANT  -  SORT  ( R4R+HH ) 

IF  (DELTA  ,LT ,  W4)  88*89 
89  INDEX- INDEX* 1 


X 
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C  COMPUTE  GRAZING  ANGLE*  PS  I*  SMOOTH-SEA  REFLECTION  COEFFICIENT,  RHO, 

C  AND  PHASE  ANGLE*  PHI*  PHASE  DIFFERENCE,  ALPHA,  ROUGH-SEA  COEFFICIENT, 
C  RUF.  AND  PATTERN-PROPAGATION  FACTOR*  FF. 

90  PSI  *  ATAN  < <H1+H2)/(6076.*R)*AK> 

CALL  SEAREF  <FMHZ *PSI « I  POL  * RHO, PH  I  > 

RATIO  «  DELTA  /  W+PHI/TUPI 

WHOLES  »  I NT  F  <  RAT  I O )  , 

RATI 01  «  RATIO  -  WHOLES 
ALPHA  *  TUP I  *  RATI 01 

RUF  r  EXP  (-8,*(PI*WH*SIN(PS1 )/W>**2) ' 

"  D  *  Dt  *  RHO  *  RUF 

FF=  SORT  ( l.+D*D+2.»D*C0S  ( ALPHA >> 

FRR  *  FF#  RFS/RSLANT 
IF  (INDEX  .GT.  2)  GO  TO  280 

C  SET  UP  4-ELEMENT  ARRAY  .OF  FIRST  TWO  POINTS  IN  INTERFERENCE  REGION  AND 
C  TWO  POINTS  IN  DIFFRACTION  REGION  (OMIT  THIS  IF  INDEX  GREATER  THAN  2), 
I  SUB* I NDEX+2 
XX  < I SUB* 1 >  *RSLANT#XR 
YY ( I  SUB ♦ 1 )»YF <  FRR) 

IF  (INDEX  *EQ«  1)  GO  TO  88 
RNM»2.#RH 

CALL  DIFFRACT  (HI  * H 2 « RNM « F MHZ . F DB  > 

FF  *10, *# ( FOB#  « 1 ) 

FRR«FF#RFS/RNM 
YY (2*1) =YF(FRR) 

XX (2* 1 )«RNM#XR 
RNM«1,1#RNM  ** 

CALL  DIFFRACT  ( H 1 «H2 , RNM*FMHZ *FDB > 

FF«10.*#(FDB*. I ) 

FRR*FF#RFS/RNM 
YY(  1.1 >*YF(FRR) 

XX ( 1 • 1) «RNM#XR 

C  DRAW  SMOOTH  CURVE  CONNECTING  THE  FOUR  POINTS*  BY  CUBIC  INTERPOLATION 
C  EQUATION. 

CALL  CURVE  (XX*  YY*  XMAX.  YMAX) 

CALL  DASHOFF 

C  RESET  VALUE  OF  FRAC  FOR  PLOTTING  INTERFERENCE  REGION. 

FRAC  ■  .9998 
GO  TO  88 

C  COMPUTE  X  VALUES  PROPORTIONAL  TO  SLANT  RANGE*  AND  Y  VALUES 
C  PROPORTIONAL  TO  DECIBELS  ABOVE  MINIMUM  DETECTABLE  SIGNAL. 

C  STOP  PLOTTING  IF  Y  GOES  ABOVE  OR  BELOW  GRID  BORDER.  RESUME  WHEN  Y 
C  RETURNS  TO  BORDER. 

280  X  «=  RSLANT  #  XR 
Y  ■  YF(FRR) 

IF ( ( Y. GT. YMAX. AND. YLAST. GT. YMAX > .OR, ( Y.LT.O. .ANO.YLAST.LT.O. ) ) 
l  30*31 

30  XLAST  -  X 
YLAST  >  Y 
GO  TO  88 

31  IF  (Y  .LE.  YMAX  .AND*  YLAST  ,L£»  YMAX)  GO  TO  43' 

IF  (V.GT.  YMAX  .AND.  YLAST  .LE.  YMAX)  40*  41 

40  CALL  INTRSECT  ( 0. « YMAX .XMAX* YMAX* XLAST * YLAST *X. Y *X I • Y I ) 

CALL  MINITAPE  (XI.YI) 

CALL  PLOT  (XJ.YI. 2) 

XLAST  >  X 
YLAST  «  Y 
M  «  3 
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GO  TO  88 

41  CALL  INTRSECT  < 0. . YMAX. XMAX. YMAX. XLAST . YLAST .X. Y ,X I ♦ Y I ) 
CALL  PLOT  (XI .YI .3) 

X2  »  XI 
Y2  »  YI 
XLAST  =  X 
YLAST  =  Y 
N1  *2 
M  a  2 
GO  TO  88 

43  IF  <  Y  .GE.  0.  .AND.  YLAST  .GE.  0.1  GO  TO  47 
IF  (Y  .LT.  0.  • AND*  YLAST  ,GE.  0.)  44*  45 

44  CALL  INTRSECT  < C. . 0. .XMAX. 0 . .XLAST . YLAST « X. Y . X! ♦ YI > 

CALL  MINI  TAPE  (XI.YI) 

CALL  PLOT  (XI .YI.2) 

XLAST  *  X 
YLAST  =  Y 
M  =  3 
GO  TO  88 

45  CALL  INTRSECT  ( 0. « 0. « XMAX. 0 .. XLAST ♦ YLAST . X * Y .XI  . YI  > 

CALL  PLOT  (XI.YI. 3) 

X2  =  XI 
Y2  »  YI 
XLAST  =  X  . 

YLAST  o  Y 


N!»2 
M  a  2 
GO  TO  88 

IF  (N1.EQ.2,  473.474 
L-X 
Y 1  =  Y  A  a  \ 

N 1  a  3 

GO  TO  88 

474  CALL  MINI  TAPE  (X.Y) 
XLAST  a  X 
YLAST  a  Y 
R  a  FRAC  *  R 
IF  'R.LT.RTEST)  880.1 
IF  (M.EQ.3)  RETURN 
CALL  MIN  I  TAPE  (X.Y) 
CALL  PLOT (X.Y. 2) 

END 


88 


880 


/ 

/ 

/ 

/ 

/ 
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/ 
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/ 
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SUBROUTINE  INVERT. ( XM IN. XMAX .NS IG « L I M.NO I . X . FT » F > 

C  THIS  SUBROUTINE  FINDS  VALUE  OF  X  THAT  RESULTS  IN  F(X>  =  FT »  BY 
C  ITERATION  BASED  ON  LINEAR  INTERPOLATION/EXTRAPOLATION  FROM  PREVIOUS 
C  TWO  TRIALS.  FUNCTION  F  MUST  BE  MONOTONIC. 

TEST  a  1 0  »•*■* ( -NS  I G ) 

FD  =  FT 

IF  (FT  .EQ.  0.)  FD  *  1. 

NO  I  =  1 

X  =  (XMAX  +  XMIN1/2. 

FI  =  FIX) 

X2  =  X 
F2  =  FI 
GO  TO  4 

1  F|  «  FIX) 

IF  (NO!  .EQ.  LIM)  RETURN 

4  TEST  1  =  ABS<  (Fl-FT  I/FO)  .  -  .. 

IF  (  TEST  1  -  TEST)  2.’  2." M" 

2  RETURN 

3  XM  =  X  - 

IF  (NO  I  .GT.  1)  GO  TO  6 

DELTA  a  (XMAX  -  XM I N > /4 . 

IF  (FT  ,LT .  FI)  DELTA  =  -  DELTA  , 

FMAX  =  F(XMAX) 

FMIN  =  F(XMIN) 

IF  (FMAX  .LT.  FMIN)  DELTA  =  -  DELTA 
X  a  X  +  DELTA 
XN  =  XM 
NOI  =  2 

GO  TO  1  . 

6  X  =  (FT  — FN)*( X  — XN ) X ( F I  — FN  >  +  XN 

NOI  =  NOI  +  1 

IF  (NOI  -  3)  24.  21  ,  24 

21  IF  ( ( ABSIF2-FT) >  -  ( ABS ( F I -F T ) ) )  23.  23.  24 

23  XN  *  X2 
FN  =  F2 

GO  TO  l  . 

24  XN  a  XM 
FN  =  F  I 
GO  TO  1 
END 


/ 


its 


,, ^ •' 
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function  essisn 

C  DEFINES  PARAMETER  S  AS  FUNCTION  OF  SI  AND  T»  SEE  RAD  LAB  VOL.13.P115. 
common/b/t.tt 
SSt  a  SI  *  SI 

ESS  a  <S1  +  (S0RT(  <  1  .-SSI)**2«-4.*SS1  *TT)-1  »  +  SSl  >/<2.0*Sl  )  >/<  1  f  +  T) 

END 


SUBROUTINE  DIFFRACT  < AHFT . THFT . RNM . FMHZ . FDB > 

COMMON/ZZX/Z1 «Z2. X.U0B1 .UDB2 

C  AHFT  IS  ANTENNA  HEIGHT  FEET*  THFT  IS  TARGET  HEIGHT.  RNM  IS  RANGE  NAUTICAL 
C  MILES.  FMHZ  IS  FREQUENCY  MEGAHERTZ.  SUBROUTINE  COMPUTES  PROPAGATION  FACTOR  IN 
C  DB  RELATIVE  TO  FREE  SPACE.  BASIS  IS  EO»  463  OF  ••PROPAGATION  OF  SHORT  RADIO 
C  WAVES. •• KERR. VOL.  13  OF  RAD.  LAB.  SERIES.  PAGE  122.  Z1  AND  Z2  ARE  ••NATURAL 
C  HEIGHTS* •  AND  X  IS  •• NATURAL  RANGE** — EQS.351  AND  358 .PAGES  96-97.  KERR. 
FACTOR  a  FMHZ  **  . 6666667/69BB. 1 03 
Z 1  <•  AHFT  «  FACTOR 
Z2  =  THFT  *  FACTOR 

X  x  RNM  *  FMHZ  *#  .3333333/102.715 
CALL  UFCN  IZ1.UDB1) 

CALL  UFCN  < Z2. UDB2) 

FDB* 1 0.99209864  +  1 0.*ALOG1 0 < X >  -  17.545497*X  +  UDB1  +  UDB2 
END 


SUBROUTINE  UFCN  <Z*  UDB ) 

C  SUBROUTINE  COMPUTES  HEIGHT-GAIN  FUNCTION  UDB  IN  DECIBELS.  BY  USING  EMPIRICAL 
C  FORMULAS  FOR  DIFFERENT  SEGMENTS  OF  FIG.  2.20.  PACE  128.  OF  • • PROPAGAT I  ON  OF 
C  SHORT  RAO  1 0  WAVES. • *KERR.  VOL.  13  OF  RAD.  LAB.  SERIES. 

C  IMPORTANT... .THIS  CURVE  IS  VALID  F Oft  HORIZONTAL  POLARIZATION  ONLY. 

IF  <Z  .LE.  .6)  I  *  2 
1  UD0»2O.#ALOG1O(Z> 

RETURN 

a  IF  (Z  .LT,  1.)  3.  4 

3  UDB  »  “4.3+BI »04#( ALOG1 0 < Z/.6 >> *# 1 *4 
RETURN 

4  UDB«19,84726*(Z*#«47  -  .9) 

END 


ij 
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autwouTiNK  curve  <xx.  yy*  xmax*  ymax> 

OIMRNSIQN  AR< A * A ) *AY ( 4 *4 ) . XX ( 4 * l ) « YY( 4 v I ) 

COMMON  /XYL/'XPR.  YPR  ' 

COMMON  /M TPE /  X{}»  V,1»*  X S  ■  Yt*  XA«  YA«  ERROR 
COMMON  AMN/  MJ *  N) 

Y(Z»  ■<  A  V(  1 4  1  )  +  E4(AY(8«t)  +  2#<  AVUill  +  Z4AYt4.ll>) 

aw  n  *»)  m  xx u . n 

AW  i  Ii3l»  \X(  1  *  J  )44£ 

API t *4 )»AR< t *3>«XX( I , 1 ) 

ARtJ*,2)»XXt2«  I  ) 

ARt£>*3)  »XX(  tl*  1)4*2 
AP<f>v4)»ARi:i,3)4XX(2.t  ) 

AR(3*8)«XX<3«U 
AM(3*3)»XX(3* 1  )4#2 
AW(  3*4  )  nAWi  3. 3)4XX<  31*  J  ) 

AN  t  I  *  I  )  ■  A»< 2. 1 >bAR(3, I) =AR( 4* 1 ) = 1 .O 
AR 1 4  *R ) »  XX<  A  «  S ) 

AR(4*3 I wXX<4* l ) 44£ 

AW (4*4)*  AR<  4*3)  4XX<  4*1) 

•— “*ATTt  » 1  >  *  YY< 

AY(2* l > »YY( 2* l ) 

AY < 3 « 1 >«YY(3*  t ) 

,  . A  Y,<  4;,^)  .  YY  <  4  .  J  )  . 

3B6  CALL  M A T ALO<  AR .  A  Y  ♦  4*1  *  0  «OT i*4  >  . .  .  . . 

XI  »  XPR  u  XMAX  \ 

Yt  «.  YPR  ■  YIXMAXJ 
-*.”*'“'"'6WL  o  (  XMAX  -  XX(3*1>)  »  .1 
IF  I  VI  .LT.  O.)  290*  291 
290  DO  292  I J  m  l,io 
XI  -  XI  -  DEL 
Y  I  »  Y  (  X  I  ) 

tF  (YJ  »GT ,  0.)  GO  TO  293 
XPR  »  XI 
YPR  «  Y I 

continue 

IF  t  Y  I  .LT,  O.)  RETURN 

CALL  (NTRSECT  < XPR . YPR , X I . Y I ,0. .0. . XMAX.O. . XO. YO> 

DEL  «■  (XO  -  XX  (3*  IV)*O.I 
XI  a  XO 
VI  .  YtXI> 

CA)  •_  PLOT  (XI  «YI  «3> 

XPR»X2=X J 
YPR«Y2oYt 
DC  300  I  I  «  I . 10 

“  XI  -  DEL  ! 

■  Y<  XI ) 

ui.eo.n  305,302 

(Yt.GT.YMAX)  303*304 

INTRSECT  (XPR*  YPR*  XJ«  YI«  0.*  YMAX *  XMAX.  YMAX,  XO.  YO> 
PLOT  (XO*  YO*  Z> 

XPR  o  XI 
YPR  »  Y I 

..return 

304  XI  «  XA  *  XPRb  XI 

. Y1  ■  YA  *  YPR  «  Yt . 

GO  TO  300 

302  IF  ( Y I  ,GT.  YMAX)  20.  21 
«*>••*•  20  CALL  INTRSECT  (XPR*  YPR* 

CALL  MINI  TAPE  (XO*  YO) 


292 


.293 


291 


X| «  YI.  0..  YMAX.  XMAX.  YMAX.  XO,  YO) 


a-"-**” 
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CALL-  PLOT  ( XO*  YO*  2) 
XPR  s  X! 

YPR  s  Y 1 
RETURN 

2!  CALL  MINITAPE  CXI.  YI > 
300  CONTINUE 
END 
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SUBPOUT INE  MATALSI  A.X.NR.NV, IDO.DET.NACT > 
DIMENSION  A(NACT«NACT) .XINACT.nIcT > 

IF  I  IDO)  1.2*1 
DO  3  I- 1 *NR 


DO  4  J«I .NR 
X< 1 , J)sO.O 
X(  I.  I  )»1.0 
NV*NR 
DETsl ,0 
NR  I »NR- I 


DO  5  K«! «NR1 
IR1 *K+i 
PIVOT  sO  ,0 
DO  6  I  sic, NR 
Z*ABSFIA< I ,K) I 
IF (Z— PIVOT)  6,6,7 
7  PIVOT *2 
IPRsI 

6  CONTINUE 

IF(PIVOT)  8.9,8 
9  DETsO.O 
RETURN 


B  IF(IPR-K)  10,11,10 

10  DO  12  JoK.NR 
Z*A< IPR, J) 

A  < IPR, J) «A(K, J) 

12  A<K,J»«Z 

DO  13  Jsl.NV 
ZsX< IPR, J) 

X< IPR, J)«X(K, J} 

13  X(K.J)*Z 
DETs-DET 

11  DET«.DET#A(K.K1 

P I VOT  « 1 . O/A ( K , K | 

DO  14  JsIRi , nr 
A<K, J)«A(K, J)«PivOT 
DO  14  IsIRi,NR 

14  A< I ,J)sA( 1, J)-A< I,K)«A(K. Jj 
DO  5  Jsl.NV 

IF(X<K.J>)  IS.  3.15 

15  X(K«J)cX(K«J) #P I VOT 
DO  16  I  «  I R 1 .,  NR 

16  X< I , J)*X( I, j)-A( I ,K)#X(K. J) 
5  CONTINUE 

I F ( A  <  NR , NR ) )  17,9,17 

17  OETsOET#A (NR.NR ) 

PIVOTs 1 , 0/A( NR, NR) 

DO  18  Jsl.NV 

X ( NR, J ) sX (NR , J) *P| VOT 

DO  18  Ksi.NRi 

I-NR-K 

SUMs 0,0 

00  19  L«I »NR1 

19  SUMsSUM+AC I ,L+l >#X<L+l , J) 

18  X< I »J)sX( l , JJ-SUM 
END 


l 
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